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ABSTRACT 


To  improve  the  stability  and  containment  time  of  a  plasma  in  a 
magnetic  field,  the  7l?sma  stabilization  for  controlled  thermonuclear 
fusion  must  be  developed.  The  methods  oi  principles  of  magnetic 
shear,  magnetic  well,  dynamic,  and  feedback  stabilization  are  of  primary 
importance  and  will  be  given  primary  consideration,  fhe  objectives  of 
this  research  are  to: 

«  review  various  stabilizing  methods, 

•  compare  the  effectiveness,  advantages  and  disadvantages  of 
some  stabilizing  methods,  and 

■  deduce  from  the  results  obtaine  l  a  new  idea  or  original 

method  for  understanding  plasma  instability  and  improving 


plasma  stabilization. 


1.  INTRODUCTION 


Plasma  stabilization  during  developing  controlled  fusion  is  one  of 
the  most  difficult  achievements  in  thermonuclear  research.  To  improve 
the  stability  and  containment  time  of  a  plasma  contained  in  a  magnetic 
field,  various  stabilizing  methods  have  been  and  are  being  developed. 

Among  them,  of  primary  importance,  are  magnetic  field  configuration 
stabilization,  dynamic  stabilization  and  feedback  stabilization.  For 
magnetic  field  configuration  stabilization,  the  plasma  may  be  made 
self-stabilizing  in  some  complex  magnetic  fields.  For  dynamic  stabili¬ 
zation,  the  high-frequency  (or  radiofrequency)  field  interacts  with 
the  entire  plasma  to  establish  a  new  dynamically  stable  equilibrium. 

For  feedback  stabilization,  the  plasma  equilibrium  remains  the  same 
and  the  time-varying  suppression  signal  interacts  only  with  the  unstable 
part  of  the  plasma.  The  reality  of  developing  a  controlled  fusion  reactor 
depends  mainly  on  the  development  of  an  appropriate  magnetic  field  con¬ 
figuration  and  the  availability  of  remote  control  and  detection  methods 
of  the  dynamic  stabilization  or  feedback  stabilization  for  a  rhermonuc lear 
plasma. 

There  are  many  methods  proposed  for  the  magnetic  field  configuration 
stabilization,  the  dynamic  stabilization  and  the  feedback  stabilization 
of  a  thermonuclear  plasma.  Only  some  of  the  main  methods  of  plasmc.  stabil 
zation  will  be  discussed,  analyzed  and  compared  in  this  research  work. 
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2.  MAGNETIC  FIELD  CONFIGURATION  STABILIZATION 

Plasma  stabilization  by  magnetic  field  configuration  consists 
mainly  of  the  following  methods:  (a)  internal  axial  magnetic  field; 

(b)  rotational  transform  devices;  (c)  screw  pinch;  (d)  magnetic  shear; 

(e)  magnetic  well-minimum  B  (open-end  system),  magnetic  well-minimum 
average  B  (clcsed-end  system);  and  (f)  multipole  configuration. 

2.1  Internal  Axial  Magnetic  Field 

To  improve  the  stability  of  a  linear  or  toroidal  pinched  plasma, 

1  2 

an  internal  axial  magnetic  field  was  early  introduced.  ’  Experiments 

show  that  the  introduction  of  an  axial  magnetic  field  in  plasma  can  delay 

instability  of  the  pinched  discharge.  Certain  ratios  of  the  axial  and 

azimuthal  magnetic  fields,  B^B^,  and  some  magnitude  of  the  pinch  ratio, 

r  /r,  (less  than  5)  must  also  be  met,  where  b  and  Bfl  are  the  axial  and 
o  z  y 

azimuthal  magnetic  fields  and  rQ  and  r  are  the  radii  of  the  discharge 

tube  and  the  pinched  plasma,  respectively. 

Interdiffusion  (or  intermixing)  of  axial  and  azimuthal  magnetic 

fields  must  be  taken  into  consideration  when  B  ard  B „  diffuse  into  each 

z  i 

other  in  the  presence  of  a  finite  plasma  resistivity,  I).  The  interdiffusion  of 
the  magnetic  fields  has  two  important  effects:,  (a)  The  redistribution  of  mag¬ 
netic  field  is  in  a  transient  state.  Aa  a  result,  the  pinched  plasma  tends  to 
become  unstable,  due  to  continuous  decrease  in  B  .  (b)  The  magnetic  field 
difru;ion  energy  loss  due  to  resistivity  (T))  in  the  current  sheath  tends  to  heat 
the  plasma,  i.e.,  the  ohmic  heating. 
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Consider  a  cylindrical  plasma  sheath  of  unit  length  and  radius  r, 
contained  by  the  axial  and  azimuthal  magnetic  fields  in  a  long  discharge 
tube.  To  evaluate  the  field  time  constant,  t,  for  diffusion  of  the  axial  field 
flux  out  of  the  cylindrical  sheath,  the  energy,  W,  of  magnetic  field  B 
(in  MKS  units)  is  given  by 


W  =  n  r 


(2.1) 


where  ^  is  the  permeability  in  vacuum.  By  differentiating  B  with  respect 
to  time,  t,  the  rate  of  the  energy  loss  is  expressed  as 


2  2 
rrr  dlT 

2p  dt 


(2.2) 


This  relation  is  equal  to  resistive  power,  P,  in  the  current  sheath  to 


heat  the  plasma 


P  =  12R  =  I21]  ~ 


(2.3) 


where  I  and  R  are  the  electrical  current  and  resistance  of  the  plasma 
and  b  is  the  sheath  thickness.  Assume  that  the  rate  of  magnetic  energy 
loss  by  diffusion  is  entirely  converted  into  the  resistive  power  of  ohmic 
heating  in  the  plasma,  from  Eqs.  (2,2)  and  (2.3),  it  follows  that 


r6  dB 
4T)gio  St 


(2.4) 


A  plasma  sheath  can  be  considered  as  a  single-turn  solenoid,  according 
to  the  Biot-Savart  law  that  current  and  magnetic  field  are  related  by 


B  =  U  I 
o 


(2.5) 
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Introduction  of  Eq.  (2.5)  into  Eq.  (2.4)  gives 

j_  agi  V! 

„2  at  '  4n 

Integrating  this  and  using  the  initial  condition  B  =  Bq  at  t  =  0,  the 

2 

time  constant  T  and  the  field  time  constant  T  related  to  B  and  B  are, 
o  ’ 

respectively 


B 


2 


B 


g-t/T 


(2.6) 


where 


T 

O 


V6 

47] 


T 


u.  r6 
o 

27] 


so  that  t  =  2t  . 

o 


It  is  obvious  that  the  diffusion  of  the  confining  magnetic  field  (or  the 
drift  of  charged  particles)  from  the  plasma  can  be  slowed  down  if  Tj  is 
small  and  r  is  large.  Tnis  result  can  be  generalized  for  various  devices: 
(a)  it  will  take  a  longer  time  for  to  diffuse  in  a  large  discharge  tube 
than  in  a  small  one,  and  (b)  the  containment  of  a  stable  discharged  plasma 
requires  rapid  initial  heating  to  attain  a  high  temperature  for  small  T|. 

2.2  Rotational  Transform  Devices 


It  is  well  known  that  the  rotational  transform  was  first  introduced 

3 

in  the  stellarator  machine.  The  stellarator  can  produce  a  rotational 
transform  because  of  its  "figure  8"  configuration.  In  the  presence  of  an 
axial  magnetic  field  with  a  rotational  transform,  the  electrostatic  fields 
and  charge  separation  due  to  the  opfositeiy  directed  drifts  of  positive 
ions  and  negative  electrons  in  the  stellarator  system  can  be  considerably 


5 


suppressed.  Magnetic  transform  and  particle  transform  can  also  serve 

4 

the  purpose  of  the  rotational  transform. 

To  compensate  for  the  ion  and  electron  drifts  in  the  curved  region 
of  c  race  track  type  of  toroidal  device,  a  corrugated  magnetic  field  can 
also  produce  the  effect  of  a  rotational  transform.  The  corrugated  mag¬ 
netic  fields  can  be  realized  by  a  suitable  spacing  of  the  winding  of 
a  number  of  solenoidal  coils  which  produce  the  axial  magnetic  field  when 
energized. 

2.3  Screw  Finch 

Screw  pinch  is  a  dynamic  pinch  wherein  an  external  axial  magnetic 
field,  B  ,  is  trapped  outside  the  plasma.  This  axial  magnetic  field  and 

b 

the  a2imuthal  magnetic  field,  Bg,  combine  to  form  a  helical  screw-type 
magnetic  field  to  confine  the  plasma. 

Early  theoretical  studies  indicated  that  the  sausage  instability 
(m  =  0)  of  a  constricted  pinch  discharge  and  the  kink  instability  (m  =  i) 
of  short  wavelength  could  be  stabilized  by  an  applied  external  axial 
magnetic  field  as  an  alternative  to  the  internal  axial  magnetic  field. 

The  possibility  of  stabilizing  a  pinched  discharge  by  an  external  axial 
magnetic  field  led  to  the  development  of  che  screw  dynamic  pinch. ^ 

The  partially  stabilized  screw  pinch  can  be  realized  by  the  appli¬ 
cation  of  the  external  axial  field  simultaneously  with  the  azi-iethal 
field  of  the  discharge.  This  *3  achieved  in  a  linear  or  toroidal  di.-.- 
charge  tube  of  the  usual  type,  modified  by  imparting  a  slight  helical 
pinch  to  the  return  current  outside  of  the  tube.  The  conductor  consists 
of  a  number  of  copper-braided  strips  (as  in  the  Tokamak  TM-3)  insulated 
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from  one  another,  arranged  along  the  length  of  the  tube  with  a  slant 
rd9/dz  equal  to  0.1-0. 2  (i.e.,  along  the  2-axis  at  radius  r  and  angle 
9).  As  a  result,  an  external  magnetic  field  is  generated  by  tha 
discharged  current,  and  the  total  magnetic  field  (both  Bg  and  con¬ 
tains  and  heats  the  plasma  in  a  helical  fashion. 

From  the  facts  that  the  region  outside  the  constricted  plasma  is  a 
partial  vacuum  and  has  appreciable  conductivity,  the  concept  of  screw 
dynamic  pinch  evolved.  To  simplify  analysis,  a  pressureleas  ionized 
gas  of  infinite  conductivity  outside  the  pinch  column  proper  is  assumed. 
The  radial  equilibrium  equation  of  magnetic  pressure  for  a  screw  dynamic 
pinch  distribution  in  the  plasma  column  is  given  by 


B 

z 


x 


d(xBQ) 

dx 


=  0 


(2.7) 


where  x  =  r/r  and  r  =  radius  of  the  linear  pinch  tube.; 
o  o  r 

Given  the  properties  of  adiabatic  invariant,  the  magnetic  fluxes, 

0 .  and  0  ,  between  a  given  point  and  the  pinch  axis  in  the  plasma  must 
b  z 

be  conserved.  Hence,  0(.  has  a  definite  functional  dependence  on  0^ . 
Then,  from  the  orthogonality  principle,  the  azimuthal  and  axial  magnetic 
fields  can  be  expressed  in  terms  of  0g  and  0^  (similar  to  the  properties 
of  the  stream  function  and  the  potential  function  in  fluid  mechanics) 


B  =  00  /r  dx 

'J  ii  o 


2 

B  =  50  /r  xdx 

z  7  0 


(2.8) 


v:hieh  satisfy  Eq.  (2.7)  of  the  pinched  plasma. 
Since 


r  dy 
o 


d0e  bK  ' 

d0  r  dx 
z  o 


(2.9) 
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Introduction  of  Eqs.  (2.8)  into  Eq.  (2.9)  yields 

B,  =  Vr  xB  (2.10) 

y  o  z 

where  V  -  d0rt/d0  which  is  a  constant  at  any  point  in  the  plasma.  This 
equation  represents  the  relationship  between  B^  and  B^  at  any  point  x 
in  the  plasma.  In  a  linear  screw  pinch,  the  return  conductor  of  the  dis 
charged  currents  is  a  helix  outside  the  plasma  tube.  If  the  radius  of 
the  pinch  tube  is  taken  as  unity,  then  the  wall  condition  of  the  screw 
pinch  becomes 

Bfl  =  V.  B  (2.11) 

9  o  z 

where 

rQx  =  1  ,  and 

v  =  (d0n/d 0  )  =  constant, 

o  y  z  , 

X=1 

Substitution  of  Eq.  (2.10)  into  Eq.  (2.7),  and  rearrangement,  give 
the  relations 


Bq  =  b  v  r  xB  , 
9  o  zo 


=  b  B 


(2.12) 


wh  ere 

» ■  o +  "l  ro>/<1  +  yl  v2)- 

at  x  =  1,  b  =  1,  and 

B  =  B 
z  zo 

on  the  wall  of  the  test  tube.  For  given  values  of  r  and  V,  the  radial 

distribution  of  B,/B  and  B  /B  at  any  point  x  of  the  screw  pinch  can 
9  zo  z  zo 


be  determined. 
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In  the  above  analysis,  the  conductivity  outside  the  pinch  column  has 
been  assumed  as  infinity.  Actually,  a  finite  resistivity  is  present.  Con¬ 
sequently,  the  interdiffusion  between  B„  and  B  as  well  as  the  decay  of 

B  will  occur  as  discussed  in  the  method  of  internal  magnetic  field 

zo 

stabilization. 

A  theoretical  study  indicated  that  a  screw  or  diffuse  pinch  with  axi- 

symmetric  configuration  operating  above  the  Kruskal-Shaf ranov  current  limit  can 

2  2 

be  stable  if  0  >  0.1  where  0  ■  2p,Qp/(Bz  +  Bg)  and  p  =  plasma  kinetic  pressure. 

For  the  purpose  of  comparison,  some  computed  results  of  a  numerical  example 

for  the  stabilized  pinch  with  external  current  and  the  screw  pincB  with 

external  current  are  shown  in  Fig.  1.  The  pinched  olasma  in  both  cases 

of  the  diffuse  pinch  is  stable.  The  stability  is  possible  for  case  (a) 

where  0  is  in  the  range  0.05  to  0.10,  while  for  case  (b) ,  0  must  be  greater 

than  0.1.  The  magnitudes  of  the  selected  magnetic  fields  B^  and  B^  are 

typical  values  and  the  pitch  is  defined  as  P  =  rB^/B^. 

Most  of  the  present  laboratory  devices  with  high  0  toroidal  p!asma 

containment  exhibit  the  flute  or  interchange  instability  (m  =  H.  Dynamic 

stabilization  of  the  m  «  1  instability  for  high  0  screw  pinch  was  studied 
8 

experimentally.  In  the  experimental  work,  two  basic  dynamic  stabilization 
methods  were  performed:  (a)  The  oscillating,  axially  uniform  high-frequency 
component  of  the  B^  field,  and  (b)  the  oscillating  axial  component  current 
methods.  For  case  (a),  special  attention  is  given  to  the  possibility  of 
working  in  resonance  with  natural  plasma  oscillation  (having  relatively 
low  frequency)  to  save  energy  input  required  by  a  high  frequency  oscillator. 
For  case  (b)  the  effects  of  dynamic  shear  on  the  flute  instability  is 
particularly  considered.  These  methods  are  similar  to  that  of  dynamic 
stabilization  to  be  discussed  later. 
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In  the  oscillating  axial  magnetic,  field  stabilization  method,  the 
oscillating  component  of  the  axial  magnetic  field  sin  a^t  can  be  pro¬ 
duced  and  superimposed  on  the  base  axial  magnetic  field  B  in  a  common 

coil  to  obtain  the  total  axial  magnetic  field,  B  =  B  +  "b  sin  a  t 

z  zo  zo  st 

where  Bzq  is  the  amplitude  of  the  oscillating  axial  field,  a*  is  the 

stabilizing  frequency  and  t  is  the  time.  The  ratio  B  / B  can  be 

zo  20 

varied  between  0.05  and  0.20.  The  range  of  u>  can  be  adjusted  from 
1-1.7  MHz.  The  growth  rate  of  the  flute  instability  is  modeiately 
suppressed  after  the  oscillating  axial  magnetic  field  has  been  applied. 

In  addition,  the  high  frequency  energy  input  could  be  saved  and  the  plasma 
column  breakdown  by  wall  contact  could  be  avoided  if  the  resonance  con¬ 
dition  between  the  oscillated  plasma  motion  and  the  natural  radial  plasma 
oscillation  could,  in  principle,  be  achieved.  In  the  oscillating  axial 
current,  the  total  axial  current,  I  =  I  +  I  sin  >  t,  is  fed  into 
the  screw  pinch  using  the  same  electrodes  for  the  high  frequency  oscillating 

I  and  the  quasi -steady  current  1  .  I  and  I  have  comparable  ampli- 

zo  zo  zo  zo 

tudes  from  15-50  kA.  The  stabilization  frequency  ou  is  about  1  MHz.  A 

reduction  of  the  m  =  1  growth  rate  is  observed  after  1  is  superimposed 

zo 

on  the  screw  pinch. 

The  growth  rate  or  instability  frequency,  f.^  =  uj ^ ^ / 2tt ,  for  the  m  -  1 
mode  instability  in  a  linear  (circles)  and  a  toroidal  (dots,  aspect  ratio 
=  1)  screw  pinch  as  a  function  of  total  axial  current  I  is  shown  in 
Fig.  2.  The  growth  rate  and  the  wavelength  of  the  m  l  instability  are' 
very  similar  in  a  linear  and  a  toroidal  high  p  screw  pinch  for  adequate', 
corresponding  plasma  parameters.  Neither  the  end  effects  nor  the  eliocts 
of  toroidal  curvature  strongly  influence  the  m  1  growth  rate,  lhe 
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Fig.  2.  Growth  rate  of  the  m  =  1  instability  in  a  linear  and  toroidal 
screw  pinch  vs  total  electric  current. 
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observed  radial  distributions  of  the  magnetic  fields  B^,  Bz  and  the  plasma 

_3 

density  n  in  a  linear  screw  pinch  at  t  =  2  ^sec,  p  =  45(10  )  Torr  and 

T  =  0.20  keV  is  given  in  Fig.  3.  It  is  seen  that  the  dynamic  effect  of 
screw  pinch  can  heat  the  ions  of  a  plasma  at  very  high  rates,  and  the 
dynamic  shear  tends  to  stabilize  the  plasma.  The  above  observed  data 

are  in  good  agreement  with  the  calculated  results  from  the  theoretical 

.  .  9,10 

analysis . 


2.4  Magnetic  Shear 


Magnetic  shear  is  appreciably  effective  in  controlling  the  gross, 

universal,  or  drift  instability  of  a  thermonuclear  plasma.  It  is  the 

electrical  short-circuit  principle  whereby  any  separation  of  charges  from 

positive  ions  and  negative  electrons  can  be  neutralized  by  short-circuiting 

11  12 

the  charges.  As  a  result,  the  plasma  can  be  appreciably  stabilized.  ’ 

The  magnetic  shear  is  produced  by  the  twist  of  magnetic  field  lines 

in  the  plasma.  The  degree  of  twist  of  the  magnetic  lines  can  be  made  to 

vary  with  distance  from  the  central  axis,  the  farther  from  the  axis, 

the  deeper  the  pitch  of  the  lines,  and  vice  versa. 

The  magnetic  shear  of  the  magnetic  field  configuration  stabilization 

can  be  produced  by  a  number  of  methods.  The  two  most  common  methods  arc 

the  reversed  magnetic  field  device  and  the  helical  magnetic  windings. 

The  influence  of  a  magnetic  shear  tends  to  make  the  radial  pressure 

gradient  positive  so  that  Suydam's  necessary  condition  of  plasma  sta- 

13  14 

bility  may  be  satisfied.  ’  In  most  case. ,  the  magnetic  shear  is  used  to 
stabilize  the  gross,  universal  or  drift  instability. The  gross  or 
universal  instability  (independent  of  plasma  configuration)  due  to  temperature 
and  density  gradients  is  synonymous  with  confinement,  and  the  drift 


13 


instability  is  the  diamagnetic  drift  caused  by  the  gradients. 

The  drift  instability  arises  from  the  Landau  damping  mechanism 
involving  particles  whose  velocities  along  the  magnetic  field  resonate 
with  wave  velocity.  Sometimes,  this  damping  mechanism  also  involves 
energy  transfer  from  the  resonant  particles  to  the  oscillating  waves. 
The  drift  instability  can  occur  under  certain  conditions;  (a)  a  large 
density  gradient  exists  in  the  plasma;  (b)  a  steep  transverse  particle 
temperature  gradient  exists;  (c)  the  temperature  gradient  opposes  the 
density  gradient  of  the  plasma;  or  (d)  a  sufficient,  local  electric 
field  exists  or  electron  current  flows  along  the  magnetic  field. 

In  the  introduction  of  the  wave  packet  model,  a  local  instability 

or  perturbation  of  the  plasma  is  assumed  to  be  composed  of  wave  packets 
$(”r,  t)  of  the  form 

$(r,  t)  =  0(r)  exp  i(E*r  -  uut)  (2.13) 


where 

$(7)  =  the  spatial  perturbation  as  a  function  of  the  space 
vector  7 , 

k  =  the  spatial  wave  number, 

uu  =  angular  frequency  of  the  plasma,  and 

t  =  Lime. 

A  general  stability  criterion  of  the  wave  packets  is  that  the  wave  energy 
of  the  wave  packets  should  not  grow  or  exponentiate  many  times  as  they  move 
through  some  region  of  local  instability. 

To  stabilize  the  drift  instability,  the  sheared  magnetic  field  in 
relation  to  the  wave  packet  has  two  main  effects:  (a)  an  effective  increase 
in  the  ion  Landau  damping  available  for  plasma  stabilization  is  through  the 


propagation  of  the  wave  packet,  or  the  expansion  of  the  potenti.il  well. 
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potential  well  can  expand  into  the  region  where  the  Landau  damping  by  the 

interaction  with  resonant  ions  is  strong;  and  (b)  the  decrease  in  wave 

energy  of  the  wave  packet,  or  the  shrinking  of  the  wave  packet  as  it 

propagates  through  some  local  region  and  gradually  vanishes. 

In  stability  analysis,  a  plasma  with  a  density  gradient  dn/dx 

(n  =  plasma  density)  in  the  x  direction,  an  unperturbed  magnetic  field 

vector  B  =  B  fcT  +  (x/L  )<T  1  and  a  Maxwellian  velocity  distribution  is 
o  L  z  s  yj 

considered,  where  B_e  =  B  =  main  axial  magnetic  field,  B  (x/L  )e 
’  0  z  z  0  s  y 

■■  ♦ 

=  shearing  magnetic  field,  Bq  =  the  scalar  magnetic  field,  e"y  and  e^  are 
the  unit  vectors  in  the  y  and  z  directions,  and  L  is  the  characteristic 
shear  or  shearing  distance. 

With  the  Boltzmann  distribution  function  fj(r,  v,  t)  =  f ^ ,  the  Vlasov 
equation  is  given  (in  Mks  units)  by 


df.  e.  ^ 

+  v-Vf,  +  — 1  (E  +  v  X  B)V.f .  =  0  , 


dt 


j  m 


j 


v  J 


(2.14' 


where  v  =  velocity  vector  in  space,  e  =  unit  electric  charge,  m  =  particle 
mass,  and  E  *  self-consistent  electric  field  strength.  The  vector  E 
satisfies  the  Poisson  equation  of  charge  (or  the  Gaussian  law) 


•  E  -  —  L  e.n 
e0  J  1  3 


in  which  en  =  permitivity  in  a  vacuum  and  n  is  evaluated  from  the  jth 


j 


species  particle  distribution  function  f ^ : 

"j  *  Jv7- 7  t)d3v- 
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If  there  is  no  initial  equilibrium  electric  field,  Eq  =  0,  the  lineari¬ 
zation  of  Eq.  (2.14)  yields  the  first-order  equation  for  f  (the  zero- 
order  equation  for  f^  is  omitted): 

df .  _  e .  e . 

-—-1  +  v  •  Vf,  .  +  -1  <v  X  Bft)V  f.  .  =  -  -L(E1  +  v  x  B.)  ■  7fn. 
at  lj  m,  0  v  lj  m.'  1  1  v  Oj 

(2.15) 

By  the  method  of  characteristics,  the  solution  of  Eq.  (2.15)  can  be 
written  as 


f.  .  - - *  (  (E.  +  v  x  B.)  •  V  f..  dt, 

lj  m  J  1  1  v  Oj 


(2.16) 


where  the  integral  is  performed  along  the  particle  orbit  in  the  lowest 
order.  It  is  assumed  that  a  time  space  dependence  of  the  perturbation 


flj  *  (fljVXP  ^ 


(2.17) 


in  which  (^  )q  is  a  constant.  By  combining  the  solution  of  the  Vlasov 
equation  and  the  well-known  Maxwell  equations,  the  dispersion  relation 


for  low-frequency  modes  of  the  plasma  can  be  given  by 

17  -1/2  A  „  ,  ....  ,2  2 


where 


7  en1/2A  T  (  {U\  kMy2 

‘-“•t 


2  bk|l  VA  1  k2  2  1 

_  - u - -  =  bkiiV. 


UU<W  -  ui  )  -  —  ~r 

e  k2  1  -  I0(b)e-” 


1  If*  h 

_J  ^ 


H'-S. 


(2.18) 


e  dn 

u  =  ~^_C  — —  *  electron  frequency, 
e  eB  ndx 


m- 


=  electron  thermal  velocity, 
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k  T.  dn 

“’i  =  =  ion  frequency. 


=  ion  thermal  velocity. 


Tg,  T  =  electron  and  ion  temperatures,  or  kinetic  temperatures, 

2  2  ^  — 
k«  =  (k  -5  k  )  =  wave  number  perpendicular  to  B  (=  B,  e  ). 

X  x  y  z  0  z  * 


k|.  =  k  +  7—  k  =  wave  number  for  shear  effects, 

II  z  K  1 

.  1.2  2  a,  12.2 

b  *  ?  ai  A3!’ 

a^  *  ion  gyroradius  or  gyromagnetic  radius, 

=  the  Alven  velocity, 

Ig(b)  =  zero  order  of  the  first  kind  of  modified  Bessel  function, 

k  ,  k  ,  k  =  wave  numbers  in  the  x,  y  and  z  directions, 
x  y  z 

In  deriving  the  dispersion  relation,  the  following  basic  assumptions 
have  been  made: 


tb:  space-time  dependence  of  perturbation  varies  as  the  behavior 
of  the  wave  packet  defined  above, 


the  characteristic  length  of  density  gradient,  r 


-1 


I-  —I 

In  dx' 


«  k^.  The  amounts  of  magnetic  shear  required  to  stabilize 
drift  instability  due  only  to  density  gradient  for  the  plasma 
/  nu  <  0  «  1  is  sought. 

the  electron  gyromagnetic  radius  is  much  smaller  than  the  wave¬ 
length. 

k,  »  k  ,  v.  «  lui/k  I  «  v 
1  z  1  1  z 1  e 

kXp  O'  1  and  w  «  0^,  where  X^  =  the  Debye  length  and 

eB 

■,  the  ion  and  electron  cyclotron  frequencies. 


Q.  -  — ,  Q 
1  m.  e 
1 


m 
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2  2 

If  T.  *  T  ,  2v./v  *  P  (=  kinetic  pressure/magnetic  pressure  of  the  plasma), 

I  c  V  A 

kl  >;>  k|j  »  kl  ”  k2»  and  c2,cj[/kjj  =  bvA  ~  veloci*;y  of  light),  then  the 
dispersion  relation  reduces  to 


(2.19) 


Based  on  the  dispersion  equations,  some  limiting  cases  for  long  and 
short  wavelengths  of  plasma  oscillation  can  be  discussed. 


2.4a  Long  Wavelength  Oscillation 


In  this  case,  b  <.<  1  and  Ig(b)e  b  —  1  -  b,  then  for  uc  =  (jj&,  Eq,  (2.19) 


becomes 


l  -  i 


_i  2  2' 

V»  .  Jlx 

2U,2 


uu(uu  -  uj  )  -  knvf  =  bkuvj  .  (2.20) 


‘||  VA  ~  ||  A 


This  equation  consists  of  three  modes  of  plasma  oscillation:  (a)  the  elec¬ 
trostatic  mode,  (b)  the  slow  Alfven  mode,  and  (c)  the  fast  Alfven  mode. 

2  2  2 

For  the  electrostatic  mode,  it  is  assumed  that  uu  «  kjjv^  and  b  s*  1 . 

The  last  term  in  Eq  .  (2.20)  is  negligible  compared  with  other  terms.  An 
obvious  solution  of  the  resulting  expression  is  uu^  =  -  as  and  ut^ 
propagate  in  opposite  directions. 

Introduction  of  +  &U)^  for  a  perturbed  frequency  uu  into  Eq.  (2.20) 

gives 


(2.21) 
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Hence,  the  perturbation  condition  of  tne  plasma  (as  («e  t  0)  is 


2b 


.2  2 

k.«  v. 


2  2  2  2  * 

1  -  2u>  /k.iv  \t) 
e  |i  A  e 


or 


_  2  ,22  ...  .22,2 
2iu  <  knv,  and  2b  .•»  ki.v./u) 
c  J|  A  ~  ||  i  e 


which  shows  that  the  unstable  electrostatic  modes  propagate  almost  perpen¬ 
dicularly  to  the  axial  magnetic  field  and  appreciable  magnetic  shear  effect 
kii  is  required  to  suppress  the  drift  instability. 


For  the  continuous  (or  transient)  Alfven  modes,  the  second  term  on 

2  2 

the  left-hand  side  of  Eq.  (2.20)  approaches  zero,  i.e.,  urfu)  -  u>e  -  kjjvA  =  0 
The  roots  of  this  quadratic  equation  represent  the  approximate  frequencies 


of  the  slow  and  fast  Alfven  modes: 

uu 


X, 


,vfi  [<•.«>*  *  *¥¥  • 


(2.22) 


Introducing  ou2  ^2  an^  ^3  ^3  ^  into  Eq. 

values  of  6ou2  and  are 


(2.20),  respectively,  the 


(2.21) 

2  2  2, 

The  perturbation  condition  for  Eqs.  (2.22)  require  for  a  Rlven 

value  of  v.,  a  relatively  small  magnetic  shear  is  needed  to  inhibit  the 

A 

drift  instability.  The  respective  frequencies  of  the  electrostatic  mode 
and  the  slow  and  fast  Alfven  modes  of  plasma  oscillation  are  represented 
by  u»l  +  6u>  (=  -  +  bu^),  uj2  F  and  uj3  +  buy 
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2.4b  Short  Wavelength  Oscillation 

In  the  oscillation  of  short  wavelength  b  >  1  and  lQ(b)e  ^  ~  (2nb) 


Then  Eq.  (2.19)  becomes 


uu 

2  +  —  |1  + 

U) 


2  r* 

L  a  **<  en  a> 

(M>>‘  +  lb  -n  - 1  HT.lt 
If*  l! 


--  0.  (2.24) 


For  the  Alfven  mode,  u/  >  4k^v2  and  the  basic  assumption  «  j uu/k | j  |  <<"  . 


e  ||  A 

Eq.  (2.24)  gives  the  unstable  Alfven  mode 


uu 


t.U  =  “ 


(2nb) 


+  ^ 


2b  k?v2 
k||VA 


(2.25) 


For  this  case,  the  effect  of  magnetic  shear  kit  oecomes  insensitive  to 

ll 

stabilize  the  drift  instability  of  the  plasma. 

2.4c  General  Shear  Stability  Criterion  for  Convective  Nodes 

To  evaluate  the  spatial  growth  of  a  perturbation  with  frequency  ur  pro 
duced  by  a  narrow  source,  the  form  of  the  relevant  approximate  solution  is 

expj  ikx^x»  w)dx 

where  k  is  a  wave  number  and  can  be  found  from  a  dispersion  relation, 
x 

This  form  is  equivalent  to  the  expression  of  the  wave  packet  (Eq.  2.11;. 

When  the  amplification  reached  by  the  perturbation  in  a  plasma  region  is 

19,20 

less  than  n  ,  the  general  stability  criterion  can  be  expressed  as 
x 

max 

-  1m  j"  k^x,  JJ)dx  <  n0  (2.2b) 

o 

in  which  exp(nQ)  represents  a  tolerable  level  of  growth  of  the  convective 


(or  transient)  modes. 
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To  illustrate  the  use  of  the  general  stability  criterion,  Eq.  (2.20) 

I 

and  the  first  of  Eqs.  (2.23)  for  the  convective  Alfven  mode  can  be  taken 
to  advantage.  As  before,  letting  u)  =  *-n  ^q*  (2*20),  it  follows 


link  = 
x 


— % 

£o  ^2 

2akll  Va 


6*. 


K1  2  2 
+  1^  L3  k||VA 


T:»tegrating  this,  and  considering  6^  as  a  parameter,  the  amplification 
of  the  growth  of  the  propagating  wave  packet  becomes 


-  1m 


max 

J 


e  w 

.  .  o  e 

k  dx  - - T~ 

x  ,.2 


K +  ^  +  Sp 


2a. kt.  v  v.  '  1^2  ^e 

i  II  e  A 


t. 

I 


d +  kti  vij  6v  t  kii  vi 


dr 


2  2  2 

for  simplicity,  assigning  the  value  2ag  ~  k^  v^,  taking  the  value  oc 
6a  corresponding  to  the  amplification,  and  performing  the  integration,  the 

I 

stability  criterion  of  the  convective  Alfven  mode  is 


L  /r<  6n  (m./em  ) 
s  otoe 


(2.27) 


where 


L  =  shearing  distance, 

s  _i 

(X  (j  n\  ^ 

— — I  =  characteristic  or  scale  length  of  density 

n  dx/ 

gradient  (or  plasma  radius),  as  defined. 


21 


2  2  2 

Recalling  the  required  condition  >  AkyV^,  the  approximation 
2 

m.  =  u)  (k.iv  /u)  )  and  carrying  out  the  integration,  the  stability 
2  e  II  A  e 

criterion  becomes 


<  2 
r  <  4 


/  v  2  /k  \2  /  m.  \  \ 

n  (  7?  |  0  ( — ~  | 

0  Vi)  \k||/  \  ome/ 


(2-28) 


It  is  important  that  with  high  g,  the  magnitude  of  magnetic  shear  needed 
for  stabilizing  the  convective  Alfven  modes  is  relatively  small. 

For  purposes  of  comparison,  for  higher  values  of  magnetic  shear  and 
larger  values  of  the  relation  exists  such  that  XjJuu/k.j  v^J  v  1,  in  the  region 
where  the  ion  Landau  damping  is  strong  enough  to  suppress  drift  instability. 
In  this  case,  JJ  —  -  ,  by  a  slight  modification  of  Eq.  (2.18),  the  stability 

criterion  for  the  electrostatic  modes  driven  by  density  gradient  is 


L  /r  <  8  (r/a,) 
s  i 


(2.29) 


which  is  corresponding  to  the  stability  criterion  obtained  by  Krall 
19 

and  Rosenbluth 

For  the  localized  region  of  the  perturbation  extending  into  the  region 

of  the  electron  Landau  damping,  Mikhailovskaya  and  Mikhailovsky  obtained, 

21 

however,  another  stability  criterion 


L  /r  <  (r/a .)  V* 

S  l 


(2.101 


which  deals  with  the  oscillations  and  stability  of  an  inhomogeneous  plasma 


contained  in  helical  magnetic  field. 
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2  2  2 

Recalling  the  required  condition  u)g  >  4k|jv^,  the  approximation 
2 

u>„  =  tu  (kiiv  /u)  )  and  carrying  out  the  integration,  the  stability 
i.  e  jj  A  e 

criterion  becomes 


t-Kfe) 


¥ 


(2.28) 


It  is  important  that  with  high  (3,  the  magnitude  of  magnetic  shear  needed 

I 

for  stabilizing  the  convective  Alfven  modes  is  relatively  small. 

For  purposes  of  comparison,  for  higher  values  of  magnetic  shear  and 
larger  values  of  the  relation  exists  such  3at  XjJuu/kjj  vJ  <  1,  in  the  region 
where  the  ion  Landau  damping  is  strong  enough  to  suppress  drift  instability. 
In  this  case,  uu  -  -  uj^,  by  a  slight  modification  of  Eq.  (2.18),  the  stability 
criterion  for  the  electrostatic  modes  driven  by  density  gradient  is 


L  /i  <  8^(r/a  ) 
s  i 


(2.29) 


which  is  corresponding  to  the  stability  criterion  obtained  by  Krall 
19 

and  Rosenbluth 

For  the  localized  region  of  the  perturbation  extending  into  the  region 

of  the  electron  Landau  damping,  Mikhailovskaya  and  Mikhailovsky  obtained, 

21 

however,  another  stability  criterion 


L Jr  <  (r/a  Kg 
s  i 


(2.30) 


which  deals  with  the  oscillations  and  stability  of  an  inhomogeneous  plasma 
contained  in  helical  magnetic  field. 
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In  the  above  wave  packet  analysis,  only  the  density  gtadient  dn/dx, 
or  dn/ndx  of  the  plasma  has  been  considered.  If  the  temperature  gradient 

dT/dx,  or  dT/Tdx  (where  T  =  T\  =  T  )  is  also  considered,  the  parameter, 

dT/Tdx  d/dx(An  T)  d(lnT).  . .  .  .  .  „  . 

dn/ndx  =  d'dx(An  n)  *  °r  d (An  n) in  t”e  X  direction  has  an  important  effect 

on  the  magnetic  shear  and  plasma  stability1^.  In  general,  for  a  large 

ion  temperature  gradient,  no  small  magnetic  shear  r/L  is  sufficient  to 

s 


stabilize  the  temperature  gradient  instability  where  r 
L  -  characteristic  shear  distance  as  defined. 


-1 


l  dn 
n  dx 


and 


<3 

From  the  wave  packet  analysis  for  both  density  gradient  and  tempera¬ 
ture  gradient  present  in  the  plasma,  the  following  analytical  results 
can  be  obtained. 

•  For  a  density  gradient  with  a  nonparallel  electron  temperature 
gradient  d  An  Tg/d  An  n  <  0,  or  a  density  gradient  with  a  small 
parallel  electron  temperature  gradient,  d  An  Te/d  An  n  <r  2,  a 

magnetic  shear  required  for  plasma  stabilization  is  r/Ls 

,  ,  ,1/3 

~  (m  /ra  ) 
e  i 

•  For  a  density  gradient  with  a  large  parallel  electron  tempera¬ 
ture  gradient,  d  An  T  /  d  An  n  >  2,  the  magnetic  shear  required 

\ 

for  plasma  stabilization  is  r/L  ~  (m  / rn.)2. 

SCI 

•  For  a  density  gradient  with  a  nonparallel  ion  temperature  gra¬ 
dient  d  An  T^/d  An  n  <  0,  or  a  large  parallel  temperature  gra¬ 
dient  d  An  T^/d  An  n  >  2,  there  is  no  magnetic  shear  related  to 
m  /m  found  for  the  plasma  stabilization. 
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2.5  Magnetic  Well 


It  has  been  experimentally  demonstrated  that  universal  or  interchange 

instability  is  considerably  inhibited  in  a  thermonuclear  plasma  contained 

22 

in  a  magnetic  well,  i.e.,  at  a  minimum  magnetic  field  —  minimum  B. 

Theoretical  analysis  also  predicted  that  interchange  instability  of  the 
plasma  can  be  contained  in  a  magnetic  well,  i.e.,  at  a  (on  the  average) 

—  23  24 

well-shaped  magnetic  field  —  minimum  average  magnetic  field  or  minimum  B.  ’ 
The  magnetic  well  may  be  conceived  as  the  flux  surface  of  magnetic  field 
(or  magnetic  surface).  In  other  words,  the  plasma  may  be  contained  in  a 
magnetic  well  where  the  field  strength  is  increased  in  every  direction 
away  from  the  plasma.  In  general,  the  magnetic  well  concerns  the  shape 
of  magnetic  field  lines. 


2.5a  Minimum 


The  concept  of  the  minimum  B  is  to  stabilize  a  thermonuclear  plasma 
in  a  magnetic  mirror  device  of  the  open-end  system.  For  the  interest  of 
a  magnetic  mirror,  if  the  requirement  of  self-consistent  equilibria  and 

plasma  potentials  in  minimum  magnetic  fields  is  satisfied,  a  stability 

criterion  for  moderate  curvature  of  axially-symmetric  magnetic  surface 
.  25 

is  given  by 


D 


(2.31) 


where  p  is  the  kinetic  pressure,  B  is  the  magnetic  field  of  the  plasma, 
and  the  integration  is  performed  along  the  element  d£  of  the  magnetic  line. 
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In  an  axially  symmetric  field,  it  is  adequate  to  take  the  magnetic 
flux  Y  for  x  as  the  radial  coordinate.  Hence  Eq.  (2.31)  becomes 


-  d  i 


1  dB  dp 
„4  dY  dY 

D 


>  0 


(2.32) 


Plasma  stability  is  feasible  if  the  pressure  is  maximum  at  the  minimum  B. 
In  other  words,  when  dp/dx  and  dB/dx  have  opposite  signs,  Eq .  (2.31)  or 
(2.32)  is  positive  and  the  stability  criterion  is  satisfied. 

From  plasma  energy  consideration,  a  stability  criterion  related  to 
the  particle  temperatures  and  the  magnetic  well  depth  for  univeral 
instability  can  be  obtained.  If  the  initial  free  energy  of  the  plasma 
is  expressed  as  the  Helmholtz  function.  A,  then  A(0)  is  given  by 


A(0)  -  (nVTe)(X2/6r)  at  t  =  0 


(2.33) 


where 


n  =  ne  =  electron  density, 

V  =  plasma  volume, 

Te  =  electron  temperature, 

X  •  wavelength,  width  of  local  instability  zone, 


r 


as  defined. 


and 


hot  the  instability  be  localized  in  a  zone  of  width  X  equal  to  a  typical 

wavelength.  A  large  number  of  charged  particles  diffuse  across  the  field 

lines.  If  the  plasma  expands  within  the  zone  \(dn/dx) ,  the  change  in 

average  energy  is  |i(dB/dx)  (where  pi  =  constant  magnetic  moment), 

oB  =  (T  +  T  ),  and  the  well  depth  VB  -  r(dB/dx) ,  then  the  initial  change 
i  e 
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in  energy,  SW,  required  for  disturbances  (in  which  plasma  expands  across 
the  field  in  every  x  direction),  is 

6W  =  IV\2ji  ~  ~  ~  nV(T.  +  T  )  (2.34) 

^  dx  dx  2  B  1  e 

r 

As  usual,  negative  6W  implies  stability.  By  combining  Eqs.  (2.33''  and 
(2.34),  the  stability  criterion  for  minimum-B  to  stabilize  universal 
instability  at  the  initial  state  is 


VB _ 1 

B  -  6(1  +  T./T  ) 
l  e 


(2.35) 


provided  that  the  oscillation  frequency  is  less  than  the  ion  cyclotron 
frequency  and  the  wavelength  is  greater  than  the  ion  gyroradius.  If 
-  T  ,  the  stability  criterion  for  the  depth  of  the  magnetic  well 
required  is  reduced  to 


VB/B  >  1/12  (2.35a) 

if  the  particle  temperatures  are  not  uniform,  for  simplicity,  let 
the  particle  temperatures  T.  ~  xdTj/dx  and  J’e  ~  xdlL,/dx  at  any  point 

x  .if  l he  thin  plasma  column,  then  the  stability  ,.  ri  t  erion  hc<  nines 


(2.36) 
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2.5b  Minimum  Average  B  (Minimum  B.  Closed-End  System) 

In  the  toroidal  system  (dlosed-end  system),  it  is  usually  found  that 
plasma  particles  diffuse  away  rapidly  across  the  magnetic  field  lines  as 
in  the  magnetic  mirror  system  (open-end  system).  To  prevent  the  universal 
or  interchange  instability  causing  plasma  particle  leakage,  the  magnetic-well 
idea  can  be  utilized  by  placing  toroidal  plasma  within  a  deep  magnetic  well. 
This  idea,  however,  cannot  be  done  by  conduction  coils  or  rods  as  in  the 
magnetic  mirror  system. 

Although  a  true  toroidal  magnetic  well  is  impossible  for  the  closed- 

end  system,  the  net  magnetic  well  or  minimum  average  B  (minimum  B)  for  the 
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stabilization  of  a  toroidal  plasma  has  been  effective.  ’  The  trend  of 
the  minimum  average  B  seems  toward  more  sophisticated  and  more  precisely 
controlled  magnetic  field  configurations. 

A  magnetic  field  with  helical  syir/try  can  provide  a  mean  or  average 
magnetic  well,  i.e.,  provide  regions  wherein  Jdi/B  decreases  away  from  a 

ii  2  2 

magnetic  axis  (or  equivalently  V  =  d  V/dY  is  negative  where  V  is  the 
volume  and  Y  is  the  flux  of  a  magnetic  surface  enclosing  the  axis  and  the 
volume).  The  integral  can  be  carried  out  along  the  element  Al  of  a  magnetic 
field  line.  In  mosc  cases,  however,  the  quantity  of  the  second  derivative  of 

the  volume  per  unit  flux,  v",  is  conveniently  evaluated  in  terms  of 


in  which  B  is  the  axial  magnetic  field  along  the  z-axis. 

Stabilization  by  negative  v"  (i.e.,  by  a  favorable  gradient  of  d//B 
or  of  average  B)  is  possible  in  toroidal  system  and  can  l>o  more  effective 
than  stabilization  by  magnetic  shear  against  the  resistive  interchange  modci. 
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Minimum  average  B  derived  from  negative  V"  of  the  toroidal  system 
against  the  universal  instability  is  about  the  same  value  as  minimum 
B  of  the  magnetic  mirror  system. 

Magnetic  wells  with  central-rod  conductor  to  carry  currents  (i.e., 

hardcore)  may  have  three  types  of  the  stagnation-point  (or  stationary-point) 

solution:  (1)  Linear  periodic-multi pole  structure,  by  using  helical  coils, 

for  instance,  1  =  2  and  4  for  shaping  field  with  normal  stagnation  points, 

and  by  superimposing  l  =  1  and  3  for  corrugating  field  to  create  favorable 

VB  regions  of  magnetic  well.  The  basic  principle  in  designing  a  negative 

V  "  system  is  that  along  any  given  magnetic  field  line  the  favorable  curvature 

should  outweigh  the  unfavorable  curvature;  ( 2 ,  Helical  flux  equilibrium, 

by  using  gross  helical  curvature  to  create  favorable  VB  regions  and  by 

using  an  axial  conductor  around  which  the  helical  equilibrium  flux  tube 

is  wound;  and  (3)  Toroidal  flux  equilibrium,  by  using  gross  toroidal 

curvature  to  create  favorable  VB  regions  of  the  magnetic  well.  For  each 

solution,  a  magnetic  scalar  potential  is  properly  assumed  to  satisfy  the 

29 

necessary  and  sufficient  conditions  of  the  problem. 

Magnetic  wells  of  the  helical  field  without  central-rod  conductor 

are  also  considered.  The  quantity  of  the  negative  V  is  evaluated  by 

27 

assuming  some  stream  Unctions  for  magnetic  flux.  From  a  practical 
viewpoint,  the  central  rod  is  inconvenient. 

2,b  Multi  pole  Configuration 

Multipole  configuration  consists  of  two  major  types  linear  multi- 
pole  and  toroidal  multipole.  The  linear  multiple  device  has  been  operated 
as  a  quadrupole,  generally  in  the  collisionless  regions,  with  plasma 
density  n  =  10^  -  10^  cm  ^  and  temperature  1-10  oV,^ 


Plasma  loss  is 


28 


mainly  across  the  magnetic  field  of  the  quadrupole.  Toroidal  multipoles 

have  been  studied  analytically  and  experimentally  from  low  p  equilibrium 

31“3r 

and  flute  instability  to  high  f3,  ballooni  »g  stability.  J  The  con- 

9  -3 

finement  of  low  density  (n  =10  cm  )  collisionless  plasma  with  T^  =  40  eV, 

T  =  10  eV  produced  by  gun  injection  or  with  T  =  1  eV,  T.  ^  1  cV  produced 
6  6  i 

by  microwave  heating  was  performed  to  determine  the  mechanisms  for  plasma 
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loss  from  a  toroidal  multipole  magnetic  field.  The  machine  has  four 
rods  carrying  current  in  one  direction  and  the  wall  cf  the  device  carrying 
the  return  current.  The  plasma  loss  was  attributed  to  various  mechanisms 
of  the  mechanical  hoops,  guard  supports,  local  walls,  plasma  fluctuation, 
diffusion,  and  the  generation  of  low  frequency  electric  fields  by  the 
guard  obstacles.  The  chief  advantage  of  the  multipole  configuration  is 
to  study  plasma  behavior  in  a  magnetohydrodynamic  (MHD)  stable  region 
or  in  the  J~d£/B  stable  region. 

To  study  the  plasma  confinement,  stability,  and  fluctuations,  the 
typical  magnetic  field  shape  in  quadrupole  and  octapole  configurations 
is  shown  in  Fig.  4.  The  magnetic  field  strength  distribution  in 
the  device  (along  the  central  vertical  line  through  the  median  plane  and 
along  the  chamber  wall)  varying  with  plasma  radius  is  given  in  Fig.  5. 

For  the  quadrupole  configuration,  the  magnetic  line  of  force  which  forms 
a  cusp,  returns  around  the  current-carrying  rods.  For  the  octapole  con¬ 
figuration,  the  four  rods  (with  magnetic  lines)  carrying  currents  are  all 
in  the  same  direction,  and  the  plasma  chamber  wall  carries  the  return 
current.  A  stabla  plasma  region  has  its  maximum  pressure  on  the'  lines 

t  =  0.  The  last  stable  or  critical  line  is  'i  .  For  stability,  the 

c 

plasma  pressure  must  decrease  from  f  =  0  to  {  =  y-  . 
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(a)  QUADRAPOLE 


(b)  OCTAPOLE 


As  pointed  out  above,  two  prin¬ 
ciples  (or  methods)  are  considered 
important  for  stabilization:  the 
principle  of  magnetic  shear  and  the 
principle  of  magnetic  well  (minimum 
B  for  the  open-end  system,  and  min¬ 
imum  B  for  the  closed-end  system 
with  short  connection  length  between 
the  favorable  curvature  and  the 
unfavorable  curvature).  In  the  toroi¬ 
dal  multiple  configurations,  the 
magnetic  shear  can  be  varied  by  super¬ 
imposing  a  toroidal  magnetic  field. 

The  depth  of  the  magnetic  well  and 
the  connection  length  can  be  altered 
by  changing  the  number  of  poles,  for 
instance,  from  octapole  to  quadruprle. 

The  superposition  of  a  toroidal 
field  produces  both  magnetic  shear 


*  =  0  STABLE  PLASMA 

*  =  CRITICAL  STABLE  LINE 

- -  =  MAXIMUM  PRESSURE  LINE 


and  finite  parallel  wave  numbers. 

In  case  the  toroidal  field  B  is  much 
smaller  than  the  average  multi  pole 


Fig.  4.  Typical  magnetic  fielo  field  B  ,  the  parallel  wave  number  k .. 

shape  in  quadrapole  and  m  H 

octapole  configuration.  is  approximately  given  by  k ..  --  (B  / B  )k 

II  f  m 

where  is  the  perpendicular  wave 
number.  The  effect  of  the  magnetic  shear  tends  to  stabilize  the  plasma  in 


a  multipole  configuration. 
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Fig.  5.  Magnetic  field  strength  distribution  in  the  device  vs  plasma 
radius . 


The  magnetic  well  depth  of  a  toroidal  multiple  configuration  can  be 

defined  as  |(Jd£ /B>Y=Y  /2  -  (JdX/B)^.  I  /  (  where  Y  is  the 

c  z '  /  z 

magnetic  flux  function  and  Y  is  the  flux  at  the  critical  line  (see  Fig.  4). 

c 

The  magnetic  well  of  the  octapole  configuration  is  usually  much  deeper 
than  that  of  the  quadrupole.  The  flux  surface  where  the  magnetic  shear 
vanishes  is  closed  to  the  Y  =  Y£  surface  when  the  aspect  ratio  of  the 
device  is  sufficiently  large. 

In  the  toroidal  multipole  devices,  drift  cyclotron  instability  and 
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interchange  (or  universal)  instability  have  been  identified.  The 
corporation  of  the  magnetic  shear  and  the  magnetic  veil  with  short  connec¬ 
tion  length  between  favorable  and  unfavorable  regions  into  the  multipole 
configuration  may  improve  the  plasma  stability. 
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3.  DYNAMIC  STABILIZATION 


The  purpose  of  dynamic  stabilization  is  to  perturb  the  equilibrium 
parameters  of  a  whole  thermonuclear  plasma  to  stabilize  it.  This  can  be 
accomplished  analytically  or  experimentally  by  assuming  (or  applying)  an 
oscillating  magnetic  field,  oscillating  electric  field,  or  microwave 

heating.  The  dynamic  stabilization  is  effective,  which  can  stabilize  many 

39  40 

modes  of  one  instability  and  other  types  of  instabilities  simultaneously.  * 
The  dynamic  stabilization  of  a  thermonuclear  plasma  may  consist  of 
the  following  proposed  methods:  (a)  Traveling  wave  and  standing  wave, 

(b)  oscillating  electromagnetic  field,  (c)  oscillating  electric  field, 

(d)  electron  temperature  modulation,  and  (e)  axial  alternating  current. 

Other  methods  of  dynamic  stabilization  may  be  developed  for  certain 
purposes  or  circumstances  once  the  basic  behavior  of  stabilization 
mechanisms  is  better  understood. 


3.1  Traveling  Wave  and  Standing  Wave 

The  dynamic  stabilization  of  linear  9  pinch  has  been  previously 
40 

proposed.  In  the  model  under  consideration  the  collisionless  plasma 
consists  of  charged  particles  which  can  be  reflected  into  it  from  its 
surface.  Every  perturbation  of  the  plasma  gives  rise  to  an  excessive 
average  magnetic  pressure  which  opposes  the  unstable  deformation  on  the 
plasma  surface.  In  addition  to  the  axial  magnetic  field,  the  applica¬ 
tion  of  a  periodic  constant  azimuthal  magnetic  field  along  the  pinch 
axis  can  produce  a  significant  effect  on  the  m  =  1  instability  of  the 
plasma. 
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The  dynamic  stabilization  of  the  0  =  1  bumpy  or  toroidal  0  pinch 

41 

has  also  been  studied.  In  the  analytical  study,  the  dynamic  stabi- 
lizing  currents  are  applied  perpendicular  to  the  axial  magnetic  field. 

(Without  the  dynamic  stabilization,  the  configuration  is  basically  un¬ 
stable.)  The  plasma  is  assumed  incompressible  because  in  the  analysis 
of  the  static  pinch  the  most  unstable  modes  of  the  plasma  are  found 
incompressible.  The  new  concept  introduced  in  the  analysis  is  to  make 
the  corrugated  magnetic  field  travel  along  the  pinch  axis  at  the  wave 

I 

velocity  vw  (where  v^  >  v^  and  v^  =  Alfven  velocity)  because  the  outward 
drift  of  the  plasma  due  to  the  curvature  of  the  torus  can  be  suppressed 
by  corrugating  the  plasma  surface.  In  other  words,  the  corrugated  mag¬ 
netic  field  and  surface  profile  are  made  to  propagate  along  the  pinch 
axis  with  a  wave  velocity  vw  by  suitably  alternating  part  of  the  current 
in  the  conduction  coils.  This  concept  is  referred  to  as  "dynamic  sta¬ 
bilization  of  the  plasma  by  traveling  wave"  (instead  of  standing  or  sta¬ 
tionary  wave  at  steady  state) .  It  may  be  assumed  that  the  frequency 
of  the  magnetic  field  should  be  greater  than  v^/\  with  the  traveling 

wave  and  less  than  v  /X  with  the  standing  wave  where  v  is  the  sonic  velo- 

s  s 

U2 

city  or  the  mean  ion  velocity  and  X  is  the  wavelength.  Dynamic  stabili¬ 
zation  by  traveling  wave  is  more  effective  than  that  by  standing  wave,  in  some 
cases.  The  idea  of  the  dynamic  stabilization  of  a  thermonuclear  plasma 
either  by  the  traveling  wave  or  by  the  standing  wave  depends  basically 
on  the  equilibrium  principle  that  all  resultant  destabilizing  forcas  at 
any  point  of  the  plasma  (contained  in  a  linear  or  toroidal  system)  must 


be  vanished. 


3.2  Oscillating  Electromagnetic  Field 


The  desirability  of  stabilizing  the  universal  instability  of  an 

inhomogeneous  plasma  in  a  magnetic  field  by  dynamic  stabilization  has 
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been  studied.  The  application  of  an  oscillating  electromagnetic 

field  in  stabilizing  plasma  is  briefly  analyzed  in  this  section.  Con¬ 
sider  a  collisionless  plasma  confined  in  the  x  direction  by  a  uniform 
magnetic  field  B^,  parallel  along  the  z  axis.  The  scale  length  of 
density  gradient  is  given  by  r  ^  =  -  dn/ndx  where  n  is  the  plasma  density 
(n  decreasing  with  increasing  x) .  The  temperature  gradient  of  the  plasma 
is  small  conpared  to  the  density  gradient.  An  oscillating  electromag¬ 
netic  field  which  consists  of  the  electric  field  E  and  the  magnetic 

field  B  can  be  expressed  in  terms  of  the  scalar  potential  0  and  the  vector 
— * 

potential  A 

■  ♦ 

— *  1  c\A  — *  — » 

E  =  -  70  -  ~  -  B  =  7  X  A  (3.1) 

where 

A  =  aQ(e1^a1(x)  +  . . .) , 

ifl  2 

0  =  aQ(e  ?>1(x)  +  ...)  +  aQ0o(x), 

0=ky+kz+Ut, 

ft  =  the  frequency  of  the  electromagnetic  field, 
c  =  the  speed  of  light, 

k  ,  k  =  wave  numbers  along  the  y  and  z  axes,  respectively, 

y  z 

0o(x),  ^(x)  =  given  functions,  and 

a  =  amplitude  of  the  electromagnetic  field  and  is  assumed 
to  be  small,  and  a^(x)  are  given  vectorial  functions. 
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The  stabilizing  effect  produced  by  the  electromagnetic  field  is  due  to 
the  drift  velocity  of  the  guiding  center  of  ions  and  electrons.  The 
dynamic  stabilization  is  effective  when  the  drift  velocity  is  of  the 
order  of  the  diamagnetic  ion  velocity  vrfi  (=  negative  diamagnefc  elec¬ 
tron  velocity  vde,  where  =  T./m.ru^.,  T.  «  T  =  ion  temperature, 
m.  =  ion  mass  and  =  ion  cyclotron  frequency). 

In  the  absence  of  the  oscillating  electromagnetic  field,  let  p 
be  the  radius  of  gyration,  v^  the  velocity  parallel  to  the  z-axis, 
and  x£  the  abscissa  of  the  guiding  center  of  a  charged  particle.  The 

motion  of  the  particle  is  not  in  resonance  with  the  electromagnetic 
field. 

In  the  presence  of  the  electromagnetic  field,  the  motion  of  the 
charged  particle  at  a  given  point  x  is  specified  as 

V  +  W  V  and  v  +  Vv  V 

where 

Ti  "V* 

t  =  (Q  +  k  v  +  k  v  )t, 

*  y  y  z.  z 

Vy,  vz  and  uJc  =  the  mean  velocities  in  the  y  and  z  directions  and  the  mean 
cyclotron  frequency  of  the  particle,  respectively, 
t  =  time,  and 

f(T1,  12)  -  the  periodic  function. 

For  given  time  t  and  particle  energy  W,  the  mean  velocities  \  the  moan 

cyclotron  frequency  of  the  particle  as  functions  of  p  ,  vn  ,  x  can  be 

c  II  c 


given  by 
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-  ,  ,  1  dW 

V  (p  ,  V,,  ,  X  )  =  -  -T— 

y  Kc’  II’  c  mu)  ox 

c  c 

~  .  .  1  dW 

v(p,v..,x)=v1,+  -  — — 
z  c  li  e  11  ra  dv|| 

-  ,  ,  2  dW 

tUvp  »  V|.,  x)  =  uu - r 

c  Kc’  II’  c  c  mu)  2 

C  3P„ 


(3.2) 


qBz 

where  uuc  =  ,  q  =  charge,  and  m  =  mass  of  the  particle. 

For  an  inhomogeneous  plasma  with  uniform  temperature  T,  the  scalar 

potential  0^  may  be  written  as 

i(k  y+k  z+u)t) 

0T  =  b(x,  0)e  y  Z 


in  which  b(x,  0)  is  a  function  of  x  and  0,  periodic  nature  in  0  with  a 
period  2tt.  Comparing  the  expressions  given  for  Eq.  (3.1),  b(x,  0)  can 
be  expanded  into  series  as 

b(x,  0)  =  bQ(x;  +  b1(x)e1^  +  b_^e  +  . , .  -  bQ(x) 

where  b^(x),  b  ^(x),  ...  are  small  compared  with  bQ(x)  and  are  of  the 
order  of  magnitude  corresponding  to  the  small  amplitude,  aQ. 

ik  x 

For  drift  modes  of  low  frequency  hi  «  hcil»  Vx>  =  h(x>e  x 
where  h(x)  is  defined  by  the  condition  Jdh/hdxj  <  l/4r^  (where 

2  h 

r.  =  [2T./m.UJ  , 1  =  mean  radius  of  gyration  of  ions),  the  dispersion 

1  11  ci 

relation  of  the  plasma  oscillations  in  the  presence  of  the  oscillating 
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electromagnetic  field  is  given  by 

L.  +  L  +  Da^  =  0  (3.3) 

i  e  ° 

where  L(L.  and  L  for  ions  and  electrons)  is  related  to  the  mean  value 
i  e 


of  the  particle 
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L  =  1 


i  .  i  \  iff  f(p  ,  v  )h^(x  )dp^dv  dx 

(U)  +  Vd)Jo(klpc)  JJJ  c  II  c  c  II  c 


u)  +  k  v  4  k  v 
y  y  z  z 


fo(pc’V||  )h2(^dp^v(|  dx£ 


where 

to  2  2  1 

■  m(a)^pc  +  Vj|  )/2TJ 


II 


.2  .2,2 
k ,  =  k  +  k 
1  x  y 


J  (k.p  )  =  the  Bessel  function  of  zero  order  of  the  first 
o  i  c 

kind,  and 

D  =  real  and  hounded  in  absolute  value  independent 
of  the  modes. 

The  plasma  is  dynamically  stabilized  if  the  real  parts  ReL^  >  0  and 

ReL£  >  0  when  |Da^|  «  1,  |u)|  <  |k  vj  ,  and  |dh/hdx|  <  1/4^.  These 

conditions  can  be  satisfied  if  the  values  of  function  v  (p  ,  v  ,  ,  x  ) 

yKc’  ||  c 

2  2  2 

for  each  set  in  the  domain,  p  <  2T/mu)  ,  v  ,  <  7/m  and  lx  |  <  r.  are 

c  c  ||  1  c'  i 

T 

spread  over  an  interval  of  the  order  of  |v^|  (  =  ^  f  =  diangnetic  velocity). 

c 

There  are  two  cases  related  to  the  dynamic  stabilization  by  the 
standing  wave  of  the  oscillating  electromagnetic  field: 

(a)  when  the  oscillating  electromagnetic  field  is  an  electromagnetic 
standing  wave  along  the  x-axis,  then 


k  =  k  =  0,  |u)  |  «  Q  »  k  (T/m  ) 

y  z  ’  c  x  e 

e 

=  the  plasma  frequency). 

Let  0^(x)  =  0,  a^(x)  =  2b  cos  k^x,  and  c 
electromagnetic  field  induced  by  an  exterral 

fa2 

0  (x)  - - ■‘r  cos  k  x. 

o  2  x 

m  c 
e 


}y  2 

2  and  U  »  (Jj  ( 'Jo  =  (4nnq  / m  ) 
P  P  e 

2  2  2  2 
k  =  ii  -  ijj  ,  for  t he 
x  p’ 

source,  it  follows  that 


1 

'? 
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Hence  Eqs.  (3.2)  are  reduced  to 


k  2~P- 

v  (p  ,  v  ,  x  )  =  — ^ — — r  J  (2k  p  )  sin  2k  x 
‘I  ‘  2m  ()2  0 

e 


y  c 


where  E  is  the  mean  square  intensity  of  the  electric  field  of  the  oscillating 
elect romegnetic  field. 

(b)  when  the  oscillating  electromagnetic  field  is  an  electrostatic 
standing  wave  along  the  x-axis:  as  before,  let  kz  =  a^(x) 

=  0Q  ( x)  =  0 

I Q/uj  .1  «  Ik  /k  1  «  1,  and  0.  (x)  =  2b  cos  k  x 
1  ci'  1  y  x1  1  x 

If  the  electromagnetic  field  is  produced  in  the  plasma,  then 
Eqs.  (3.2)  become 


v  (p  ,  v  ,  x  )  = 
y  Kc  il  *  c 


2  2 

2q  E  k 

- r — ^  J^(k  p  )  cos  2k  x  . 

2n  o'  xrc  x  c 

mu)  Q 
c 


3.3  Oscillating  Electric  Field 


Oscillating  high-frequency  electric  field  applied  for  dynamic  sta¬ 
bilization  of  an  inhomogeneous  plasma  consists  of  two  cases: 

(a)  Uniform  high-frequency  electric  field,  when  the  skin  laver 

6^  of  the  field  is  equal  or  greater  than  the  transverse  dimensions  of 
the  plasma  L^,  i.e.,  when  6^  >  Lj^,  the  electric  field  is  assumed 
uniform. 

(b)  Non-uniform  high-frequency  electric  field,  when  the  skin  layer 
depth  <  L|,  the  electric  field  is  considered  non-uniform  in 
the  plasma. 

In  both  cases  the  stabilizing  effects  depend  on  the  plasma  parameters 
and  the  frequency  and  strength  of  the  electric  field. 
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A  thermonuclear  plasma  confined  by  a  magnetic  field  in  a  closed-end 

device  is  subject  to  both  hydroraagnetic  and  drift  instabilities.  These 

instabilities  can  result  in  anomalous  diffusion  of  the  plasma  across  the 

magnetic  field.  High-frequency  dynamic  stabilization  of  hydromagnetic 

40-45 

instability  has  been  studied;  dynamic  stabilization  of  drift  (gross 

or  universal)  instability  by  oscillating  electric  field  has  also  been 

.  46-48 
considered. 


3.3a  Uniform  High-Frequency  Electric  Field 

The  stabilization  of  collisionless  drift  instability  in  an  inhomo¬ 
geneous  plasma  is  feasible  by  the  application  of  a  uniform  high-frequency 
electric  field  E(t)  in  the  direction  parallel  to  the  magnetic  field, 
which  may  be  called  the  "E-stabilization": 

E(t)  =  Eq  cos  Qt.  (3.4) 

where  E(0)  =  Eq,  t  is  the  time  and  (2  is  the  frequency  of  the  applied 
electric  field.  When  a  high-frequency  electric  field  of  sufficient 
amplitude  operates  in  the  plasma,  the  oscillations  of  electrons  in  the 
field  become  important.  If  a  drift  instability  is  excited  in  the  plasma, 
the  plasma  electrons  are  subject  to  an  additional  force  in  ttie  direction 
parallel  to  the  magnetic  field  due  to  the  presence  of  the  high-frequency 
electric  field,  E(t)  and  E^(t^,  x).  This  force  is  equal  to  the  electron 
charge  times  the  mean  values  of  the  electric  field,  e(E  +  E^),  where  E^ 
is  the  part  of  high-frequency  induced  electric  field.  The  parameter 
d Un  T)/d ( in  n)plays  a  very  important  role  in  the  determination  of  sta¬ 
bilizing  effect  for  different  frequency  regions  related  to  the  parameter 
of  ion  speed  ratio,  (here  T  =  plasma  temperature  and  n  =  plasma  density 


as  defined). 
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When  this  additional  force  is  in  phase  with  the  force  due  to  kinetic 

pressure  of  the  plasma,  the  electric  field  in  the  drift  wave  is  increased 

and  so  is  the  frequency,  U).  Under  the  conditions  (T^/nu)2  «  w/k^  «  (T^/nO2 

k  2  2 

and  L  =  (k, /uu,  )  (2T . /m.  )  of  the  intermediate-frequency  case  (k.  =  k 
si  it  11  x  x 

2 

+  ky,  =  Larmor  frequency),  as  uu  is  increased,  the  Landau  damping 

by  the  electrons  is  intensified.  Therefore  the  collisionless  drift 

instability  of  the  plasma  tends  to  be  stabilized.  The  range  of  values 

of  the  parameter  d(X,nT)/dUn  n)  in  which  high-frequency  stabilization 

occurs  is  maximal  for  ->  0  and  approaches  zero  when  -  ».  For  the 

low-freqeuncy  case,  uu  «  k^T^rn^  the  instability  region  is  the 

same  as  in  the  absence  of  high-frequency  electric  field  d  (£n  T)/d(£n  n) 

>  2(1  +  6)  where  6  =  ^[l  -  I^cJ/2) /Iq(cJ/2)]  and  lo(pJ/2),  l^pJ/2) 

are  the  modified  Bessel  functions. 

For  the  high-frequency  case,  uu  »  kz(Te/me>  ,  if  the  conditions 
l  r 

j  <•'"  k  (  (T  /m  K  and  k  «  k  t.(m  /m. )  are  giver,  and  the  requirement  of  l  lu- 
oDiii  z  oiei 

parameter  d(-»in  T) /d(£n  n)  <  (l/£  -  D  is  satisfied,  the  dynamic  stabilization 

of  a  collisionless  inhomogeneous  plasma  by  an  oscillating  uniform  high-frequency 

electric  field  with  modest  amplitude  can  be  achieved.  The  value  of  kQ  is 

defined  as  k  =  (n  T)-1d(n  T)/dx,  where  n  and  T  are  the  plasma  density  and 
o  o  o  o 

kinetic  temperature  at  equilibrium  state. 

The  electron-ion  collisions  in  a  collisional,  inhomogeneous  plasma 

resulting  in  the  excitation  of  drift  instability  is  known  as  the  drift- 

49 

dissipative  instability.  When  a  dispersion  relation  is  derived  on  the 
basis  of  the  usual  hydrodynamic  equations,  the  collision  term  m^v^y 
due  to  electron-ion  collisions  must  be  taken  into  account  in  the  equation 
of  electron  motion  parallel  to  (where  m^  is  the  electron  mass,  v^  is 
the  electron  velocity  and  v  is  the  electron-ion  collision  frequency  in 
the  z  direction)..  It  can  be  shown  from  solutions  of  the  dispersion 
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equation  that  the  drift-dissipative  instability  can  also  be  stabilized 

by  a  high-frequency  electric  field.  The  amplitude  of  the  uniform  high- 

frequency  electric  field  required,  however,  is  higher  than  that  for 

47 

the  collisionless  drift  instability  obtained. 

3.3b  Non-Uniform  High-Frequency  Electric  Field 

As  discussed,  if  the  skin  layer  depth  6^  of  the  high-frequency 

electric  field  is  comparable  with  the  transverse  dimension  LA  or  plasma 

radius,  the  high  frequency  electric  field  can  be  considered  as  uniform 

47 

and  the  effect  of  the  magnetic  field  can  be  neglected  If  oA  «-  LA 
the  high-frequency  field  skin  effect,  however,  must  be  taken  into  account 
The  high  frequency  electric  field  is  then  non-uniform  and  the  effect 
of  the  magnetic  field  component  becomes  important.  In  such  case,  the 
stabilizing  effect  of  the  non-uniform  high-frequency  electric  field 

48 

gives  rise  to  two  additional  independent  mechanisms  of  stabilization 
(a)  the  gravitational  mechanism  produced  by  the  averaged  high-frequency 
field  forces  in  an  unperturbed  plasma  state,  (b)  the  oscillating 
magnetic-field  mechanism  of  the  time-dependent  component  of  the  wave 
field  vector  projected  in  the  axial  magnetic  field  B  direction. 

The  gravitational  mechanism  of  stabilization  may  be  called  the  "f-stabili- 
zation,"  and  the  oscillating  magnetic  field  mechanism  of  stabilization 
may  be  called  the  "B-stabilization . "  Thus,  the  effects  of  the  B-stabili- 
zation,  E-stabilization  and  g-stabilization  of  drift  waves  by  a  non-uniform 
high-frequency  electric  field  can  occur  concurrently  on  a  weekly  inhomo¬ 
geneous  plasma  confined  in  a  constant,  main  magnetic  field,  B  . 
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Consider  a  collisionless,  inhomogeneous  plasma  confined  in  a  strong 
constant  magnetic  field  along  the  z-axis.  The  applied  non-uniform 
high-frequency  electric  field  can  resolve  into  two  components: 


E  =  E  (x)  sin  ftt,  and 
z  o 

c  dEo(x) 

B  =  -  -x  — r -  cos  Qt 

z  11  ox 


(3.5) 


where  Eq(x)  is  the  amplitude  of  the  applied  electric  field  E,  B^  is  the 
induced  magnetic  field  component  parallel  to  the  main  magnetic  field 
,  and  Q,  t  and  c  are  as  defined  before.  The  plasma  density  n  and 
temperature  T  vary  with  x  radially. 

The  averaged  high-frequency  field  force  which  influences  the 
g-stabilization  in  an  unperturbed  plasma  state  is  given  by 

F«  ■  -A  -  -  Va  ■  -  V  ^2^  <3-*> 

a 

where 


g  =  the  gravitational  acceleration, 

Y  =  the  high-frequency  potential  function, 

E  =  E  (x)  in  the  z  direction,  and 
zo  o 

a.  =  e,  i  for  electrons  and  ions  in  the  plasma  respectively. 

The  high-frequency  field  forces  act  mainly  with  the  negative  electrons. 

In  the  steady  state,  however,  there  is  charge  separation  together  with 

a  potential  electric  field  E  =  -  V0  (see  Eqs.  3.1),  or  =  -  o0/dx. 

As  a  result,  the  positive  ions  are  also  contained. 

The  conditions  of  a  collisionless,  inhomogeneous  plasma  require 

,  -1  1  dn 

that  the  dimensions  of  particle  density  gradient,  r  =  -  ~  and 


42 


■*1  1  dT 

particle  temperature  gradient,  (r  *)  =  -  —  — ,  and  the  skin  layer  depth  6^  of 

the  high-frequency  electric  field,  are  large  compared  to  the  Larmor  radius 
of  the  ions. 

From  the  oscillating  electric  and  magnetic  fields  and  the  field  force 
produced  by  the  non-uniform  high-frequency  electric  field  and  the  con¬ 
ditions  required  for  a  coilisionless ,  inhomogeneous  plasma  confined  in 
a  constant  magnetic  field,  the  dispersion  relation  can  be  derived  and 
solutions  of  the  dispersion  equation  can  be  obtained.  In  the  absence  of  a 
high-frequency  electric  field,  the  plasma  is  stable  only  against  a  fast 
ion  acoustic  wave  when  the  parameter  dUnT)/d(£nn)  >  0.  The  plasma  is 
unstable  for  d(£nT)/d(£nn)  <  0.  When  the  high-frequency  electric  field 
is  applied,  the  stable  region  broadens  toward  negative  values  of 
d(£nT)/dUn  n).  The  effect  of  stabilization  is  the  increased  frequency 
Oi  which  accompanies  the  applied  high-frequency  electric  field,  as  a 
result,  the  Landau  electron  damping  is  intensified.  Thus,  the  stabilizing 

effect  produced  by  the  high-frequoncy  field  pressure  i?  associated  with  the 

8a 

Doppler  shift  due  to  drift  velocity  v  =  -*  (where  J)  =  qB  /m  )  of  the 

Cl  U)  ( x  Z 

a 

charged  particles,  i.e.,  the  mechanism  of  the  g-stabilization. 

In  general,  the  stability  of  a  collisionless,  inhomogeneous  plasma 
stabilized  by  a  non-uniform  high-frequency  electric  field  depends  on 
the  plasma  parameters  and  the  frequency  and  strength  of  the  electric  field. 

An  effective  stabilization  of  the  plasma  by  a  non-uniform  electric  field 
may  be  achieved  near  the  resonance  condition. 

Analytical  work  to  predict  the  effects  of  dynamic  stabilization  on 
drift  instability  of  an  inhomogeneous  plasma  by  an  oscillating  electric 
field  parallel  to  the  magnetic  field  h3s  been  discussed  above.  Very 
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recently,  experimental  work  to  observe  the  effect  of  dynamic  stabili¬ 
zation  of  the  plasma  by  the  oscillating  electric  field  has  also  been 
reported. in  the  experimental  work  reported,  dischage-type  helium 
plasma  and  Q-;nachine  cesium  plasma  have  been  used.  Plasma  diagnostics 
are  made  by  means  of  radially  movable  electrostatic  (or  Langmuir),  and 
magnetic  probes.  Electron  density,  electron  temperature,  magnetic  field 
strength,  modulation  frequency,  wave  amplitude,  ac  voltage,  etc.,  are 
measured.  The  decreases  of  drift  wave  amplitude  and  wave  spectrum  in 
several  modes  (m=l,m=2,  ...)  prove  the  stabilizing  effect  of  the 
electric  field  at  relatively  low  and  high  frequency  regions. 

A  proposed  schematic  diagram  for  dynamic  stabilization  of  a  dis¬ 
charged  plasma  by  a  high-frequency  electric  field  is  shown  in  Fig.  6. 
Basically,  the  system  consists  of  two  cathodes  (with  a  tubular  anode 
between  them)  in  a  magnetic  field  parallel  to  the  axis  of  the  anode 
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Fig.  6.  Schematic  diagram  for  dynamic  stabilization  of  a  d ischarged plasma 
by  high-frequency  electric  field. 
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tube.  Deuterium  gas  is  discharged  through  perforated  tantalum  cathode 
at  the  right,  and  plasma  emerges  from  a  tabular  cathode  at  the  left. 

Ions  produced  in  the  discharge  are  accelerated  back  and  forth  between 
the  two  cathodes,  where  the  ions  produce  secondary  electrons.  These 
electrons  are  prevented  by  the  magnetic  field  from  reaching  the  anode, 
and  they  in  turn  produce  more  ions  in  the  plasma.  Neutral  particles  are  removed 
by  the  divertor.  A  sinusoidal  ac  potential  with  variable  amplitude  and 
frequency  is  applied  near  the  left  end  of  the  plasma  column.  The  average 
current  of  the  plasma  is  measured  by  ammeter  A,  Fig.  6.  Plasma  diagnostics 
can  be  made  by  means  of  movable  electrostatic  probe,  magnetic  probe, 
spectroscopic  measurement,  visual  analyzer  and  microwave  techniques. 

3.4  Electron  Temperature  Modulation 

The  effect  of  a  small,  periodic,  time-dependent  modulation  of  the 
electron  temperature  on  the  drift  waves  of  a  collisionless,  inhomogeneous 
plasma  is  considered.  The  energy  source  required  for  electron  temperature 
modulation  can  be  a  periodic,  intensity-varying  microwave  to  heat  the 
plasma.  If  the  transit  time  of  a  particle  through  the  plasma  device  is 
short  compared  to  the  modulation  period  of  microwave  heating,  the  elec¬ 
tron  temperature  will  follow  the  modulation  of  the  energy  source.  Energy 
transfer  from  the  electrons  to  the  ions  in  the  plasma  will  be  small  if 
the  electron-ion  equipartition  time  is  long  compared  to  the  modulation 
period.  Ohterwise,  the  energy  transfer  from  the  electrons  to  the  ions 
will  be  large. 
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If  the  plasma  has  a  Maxwellian  velocity  distribution,  the  moduiati 
electron  temperature,  Tg(x,  t),  may  be  expressed  as 


T  (x,  t)  =  T  (1  +  a  cos  ft  t) 
e  eo  m 


(3.7) 


where  T  is  the  peak  electron  temperature,  a  is  the  parameter  a  «  1, 
eo 

and  ft  is  the  modulation  frequency  of  the  energy  source  such  as  micro- 
m 

wave  heating.  By  using  the  Vlasov  equation  of  electron  motion  and 

giving  a  small  perturbation  of  the  equilibrium  parameter,  a  dispersion 

52 

relation  can  be  derived.  For  the  plasma  with  weak  density  gradient, 
temperature  gradient  and/or  axial  current  in  the  limit  of  -.anishing  Larmor 

radius,  the  stability  criterion  for  different  modes  of  drift  wave  can  be 

,  52,53 
given  by 


where 


2  2tt(uT)2  f*|k  | u  /d  In  T 


.  2 

k  v  v. 
z  e  l 


d  In  n 


dinT 
_ ec 

d  jtn  n 


,  d  In  T. 

_1  _ _ l 

2  d  In  n 


(3.8) 


,  k  v 
*  _  y  e  2.  dn 

^e  2 |ft  f  n  dx  ’ 

r  p 


v  =(2T  /m  ) 2 , 
e  eo  e  ’ 


v.  =  (2  T./m.)2, 
l  i  i 


ft  =  qB/m  , 
ce  e 


u  =  macroscopic  electron  velocity,  u  «  v^, 
k  ,  k  =  the  wave  numbers  in  the  v  and  z  directions,  and 

y  ^ 


T  ,  n(n^  =  ni),  etc.,  =  previously  defined. 
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Equation  (3.8)  is  valid  as  long  as  a  «  1  and  v^  «  j— j  «  v  .  It  is 
seen  that  the  plasma  parameters,  d(£n  Tg0)/d(jtn  n)and  dOfcnTp/dffn  n) 
remain  to  play  the  important  role  to  stabilize  the  drift  waves  by  electron 
temperature  modulation. 

The  theoretical  results  show  that  the  dynamic  stabilization  of 
drift  waves  by  electron  temperature  modulation  for  a  collisionless 
inhomogeneous  plasma  is  possible  if  the  stability  criterion  E^.  (3.8) 

ic 

is  satisfied  and  one  of  the  beat  frequency,  uj  =  +  it  -  di  is  of  the 

b+  —  m  e 

order  of  | k_^ | where  the  ion  Landau  damping  becomes  intensified. 


3.5  Axial  Alternating  Current 


Dynamic  stabilization  of  collision-dominated  drift  waves  by  means 
of  an  axial  alternating  current  (or  ac  current)  is  studied.  With  the 
collision-dominated  region,  an  inhomogeneous  plasma  has  a  uniform  mag¬ 
netic  field  B(=  B  )  in  the  z  direction  and  the  density  n  ana  the  temperature 
z 

T  varying  in  the  x  direction.  When  an  axial  electric  voltage  is  applied 

to  the  grid  of  the  plasma  tube  in  the  presence  of  an  ac  current,  the 

drift  waves  of  the  plasma  at  low  frequencies  can  be  suppressed. 

Drift  waves  in  a  cesium  plasma  are  stabilized  experimentally  by 

applying  an  external,  ac  current  parallel  to  magnetic  fields  through 
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a  meshed  grid  immersed  in  the  plasma.  The  results  of  the  dynamic 
stabilization  shown  in  the  experimental  data  may  be  interpreted  by  the 
facts:  (a)  The  bias  effect  due  to  rectifying  characteristics  of  a 

plasma  sheath  when  ac  voltage  or  current  is  applied  to  the  grid, 
causes  the  electron  current  to  increase  because  of  rectifying  charac¬ 
teristics  of  the  electron-rich  plasma  sheath,  and  (b)  some  non-linear 
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Landau  damping  effect  due  to  couplings  between  drift  waves  and  external 
perturbations  of  the  applied  current  when  the  applied  current  frequency 
Q  is  closed  to  drift  wave  frequency,  and  strong  coupling  between  them 
is  observed. 

The  simple,  analytical  result  for  dynamic  stabilization  of  collision- 
dominated  drift  waves  in  the  presence  of  an  ac  current  is  also  obtained 
when  the  applied  current  frequency  Q  is  comparable  with  drift  wave 
frequency. ^  A  macroscopic  electron  velocity  uq(x,  t)  parallel  to  the 
magnetic  field  can  be  produced  by  the  ac  current  flowing  in  the 
plasma.  From  the  equation  of  electron  motion,  the  linearized  equation 
of  continuity,  the  energy  balance  equation  of  the  system,  and 
simplification,  the  stability  condition  for  low-frequency  drift  waves  is 
given  by 


where 


3u2  d  In  T 
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2  d  £n  n 
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d  In  T 
_ eo 
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(3.9) 


.  1,  2  2  1,,  2  .  ,2,  2 

b  =  — k I  a .  =  — (k  +  k  ) a . , 
2  1  l  2  x  y  i 


=  the  ion  gyramagnetic  radius, 


v  =  (2  T  /m  )\ 
e  e  e 

V  =  the  electron-ion  collision  frequency,  V  -*  ■  |k  |v  and 

o  =  the  initial  or  unperturbed  state  of  the  plasma.; 

It  is  seen  from  Eq .  (3.9)  that  the  combined  effects  of  finite  ion 
Larmor  radius,  the  temperature  gradient  and  density  gradient  of  ions  and 
electrons,  the  currer-  gradient,  and  the  electron  inertial  influence  are 
all  involved  in  the  stability  condition  of  collision-dominated  drift  waves,. 
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4.  FEEDBACK  STABILIZATION 


Plasma  stabilization  by  feedback  methods  has  been  introduced  quite 
56-68 

recently.  The  selected  and  proposed  feedback  methods  to  stabilize  a 

thermonuclear  plasma  are:  (a)  electrostatic  feedback;  (b)  neutral  atom 
beam  feedback;  (c)  modulated  microwave  energy  sources,  and  (d)  modulated 
uniformly  distributed  electron  sources.  Some  new  concepts  of  oscillating 
or  helical  magnetic  fields  for  feedback  stabilization  of  a  thermonuclear 
plasma  can  be  developed  similarly. 


4.1  Electrostatic  Feedback 


Feedback  stabilization  consists  of  sensing  the  plasma  instability  at 
any  point  in  the  plasma  by  sensor  devices  and  then  feeding  back  an 
accurate  signal  of  proper  amplitude  and  phase  in  order  to  suppress  the 
instability.  A  number  of  feedback  stabilization  experiments  have  been 
conducted  with  electrostatic  (or  Langmuir)  probes  as  sensors  as  well  as 
control  elements.  Such  experiments  have  proved  successful  in  suppressing 
collisional  drift  instability  of  the  plasma.  The  theory  of  electrostati- 
feedback  stabilization  of  drift  waves  to  interpret  the  experimental 


results  has  also  been  advanced 
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Various  electrostatic  feedback  t  chniques  have  been  iluve Loped  for 

plasma  stabilization.  By  modulating  a  parallel-electron-current  sink, 

for  example, (i.e. ,  modulating  and  drawing  electron  current  from  an 

electrostatic  probe  parallel  to  the  axial  magnetic  field)  the  collisional 

drift  wave  is  stabilized  with  negative  feedback  in  the  unstable  region 

of  the  potassium  plasma  at  T  =  0.23  ev,  n  =  (10^)cm  \  dn  /n  dx  =  -  1.25  cm  * 

o  o  o 

and  the  critical  magnetic  field  =  2.19  kilogauss  (kG).  The  drift 
instability  is  the  m  =  2  azimuthal  mode.  The  experimental  results  obtained 
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for  constant  feedback  gain  varying  with  phase  (angle)  shift  of  the 
applied  probe  (suppressor)  voltage  for  various  magnetic  fields  are  shown 
in  Fig.  7(a-c).  Figure  7(a)  shows  the  plasma  density  or  instability 
amplitude,  Fig.  7(b)  gives  the  change  of  plasma  density  at  the  center 
of  the  plasma  column,  and  Fig.  7(c'  vhows  the  instability  frequency  as 
functions  of  feedback  phase  shift  in  both  stable  and  unstable  regions 
respectively.  The  theoretical  instability  growth  rate  and  frequency 
of  the  plasma  are  also  shown  in  the  curves  of  Fig.  7(a)  and  (c).  The 
feedback  gain  varying  with  the  plasma  density,  characteristic  length  of 
the  density  gradient  and  instability  frequency  of  the  plasma  is  shown 
in  Fig.  8(a-c).  The  radial  distribution  of  equilibrium  plasma  density 
n^,  plasma  density  or  instability  amplitude  n,  and  required  minimum 
probe  (suppressor)  voltage  V  for  feedback  stabilization  is  given  in 
Fig.  9.  The  plasma  density  at  the  center  of  the  plasma  column  as  a 
function  of  magnetic  field  strength  with  and  without  feedback  stabilization 
of  collisional  drift  waves  is  shown  in  Fig.  10. 

With  a  suitable  feedback  circuit  arrangement,  stabilization  of  the 

m  =  2  aximuthal  mode  can  be  achieved  up  to  B  =  1.2B  (B  =  critical  magnetic 

c  c 

field  strength).  Thus,  the  electrostatic  feedback  (negative  or  positive) 
demonstrates  that  collisional  drift  instability  driven  by  density  gradient 
can  be  controlled  from  penetration  of  feedback  signa's  into  tne  plasma 
interior  with  phase  shift. 

4.2  Neutral  Atom  Beam  Feedback 

The  application  of  feedback  stabilization  by  means  of  electrostatic 
P'-obe  or  material  suppressor  is  limited  to  plasmas  of  moderate  density 
,..nd  temperature  because  an  electrode  cannot  resist  very  high  temperature 


Constant  feedback  gain  vs  phase  shift  of  the  applied  probe 
(suppressor)  voltage  for  various  magnetic  fields-  (a)  instability 
amplitude;  (b)  change  of  density  at  center  of  plasma  column;  and 
(c)  instability  frequency  after  subtraction  of  Doppler  shift  due 
to  radial  electric  field. 


CHARACTERISTIC  LENGTH 


Feedback  gain  vs:,  (a)  plasma  density;'  (b)  characteristic  length 
of  the  density  gradient;  and  (c)  instability  frequency. 
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Fig.  9.  Radial  distribution  of  plasma  density,  nG,  instability  amplitude 
(density),  n,  and  required  minimum  probe  (suppressor,  rms),  vol¬ 
tage,  V,  for  feedback  stabilization. 


MAGNETIC  FIELD  STRENGTH,  B,  kG 

Fig.  10.  Plasma  density  at  the  center  of  the  plasma  column  as  a  function 
of  magnetic  field  B  with  and  without  feedback  stabi lization  of 
collisional  drift  waves. 
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of  the  plasma.  For  a  thermonuclear  plasma,  the  injection  of  a  neutral  atom 
beam  to  provide  feedback-controlled  volume  sources  of  particles  and  momen¬ 
tum  density  is  a  basically  different  approach  from  the  electrostatic  feed¬ 
back.  Any  perturbation  imparted  to  the  plasma  by  a  neutral  atom  beam  can 
be  sensed  by  optical  or  microwave  techniques. 

By  using  the  equations  of  motion  of  ions  and  electrons  (i.e.,  two-fluid 
equations  for  the  plasma),  and  the  equation  of  continuity  with  source  terms 
of  density  and  momentum,  after  linearization  and  simplification,  a  quadratic 
equation  of  the  dispersion  relation  for  frequency  uu  is  given  by^ 

k.v  v 

o(uu  -  uu.  -  uur)  +  ioj  (uu  -  o  -  iuur  +  +  -^0,  =  0  (4.1) 

if  s  e  fuufrf 

c 

where 


uj  =  k.v,  , 
e  1  de’ 

uu.  =  k.v,. , 
l  1  di’ 


o  = 
s 


'kz\2 


Mn  T\* 
o 


o 


M  ’ 

T. 

vdi  =  =  ion  diamagnetic  drift  velocity, 


Vde  =  qBr  =  e^ectron  diamagnetic  drift  velocity, 
o^  =  S/n^  =  feedback  frequency  or  gain, 


v  =  neutral  atom  beam  velocity, 
x 


-1 


M  =  neutral  atom  mass, 
dn 


n  dx’ 
o 
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=  equilibrium  plasma  density, 

=  neutral  atom  density, 

S  =  neutral  atom  beam  source, 

T)  =  plasma  resistivity  related  to  uu 

s 

T\  ,Tg»q ,B,k  ,kz,etc.,  have  been  defined. 

For  it  »  | uu  |  and  uu  =0,  Eq.  (4.1)  reduces  to 

S  G  t 

2 

uu  -  uu(uu.  -  im  )  -  iuu  u)  =0. 
i  s  s  e 

This  has  the  two  solutions 

it  =  it  +  iit  (it  -  it  )/tt  ,  and  uu  =  uu.  -  it  -  it  (4.2) 

ae  eeis  bie  s 

The  purpose  of  the  analysis  is  to  use  the  neutral  atom  beam  to  stabilize 

xa  without  destabilizing  tt^. 

There  are  three  basically  different  stabilization  mechanisms  which 

are  characterized  by  the  phase  (angle)  shift  9  and  the  relationships  between 

v  and  rt  (or  normalized  beam  momentum  v  / rit  )  of  the  neutral  atom  beam 
x  e  x  e 

feedback  system: 

•  Plasma  density  smoothing,  i.e.,  smoothing  density  perturbation 
with  the  neutral-beam  injected  plasma.  For  this  case 

)v^|  «  J ruje I  anc*  9  -  180°.  The  marginal  stability  requirement 
in  Eqs.  (4.1)  and  (4.2)  is  if)..  =  uu  (uu  -  a)  )/ix.  For  x,  »  I  x  I 

t  S  Q  S  '  O.  ' 

»  Iw  I,  this  becomes  uu.  =  uu  (uu  -  uu.)/ju  =  -  v  =  the  instability 
1  f  f  e  e  i  s 

growth  rate. 

•  Equivalent  (or  simulating)  minimum  B  configuration  |v  |  »  |  r j  | 

X  G 

and  h  =  0°.  The  marginal  stability  condition  is  =  -  j  (i  -  i,),r/v  . 

f  c  e  l  x 

The  pressure  exerted  by  the  momentum  of  the  neutral  beam  for  this 
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case  becomes  important  and  produces  an  equivalent  or  simulating 
minimum  B  effect. 

•  De-energization  of  resistive  drift  waves.  For  this  case 

|v^j  »  0  =  -  90°.  The  marginal  stability  requirement 

is  the  feedback  gain  Ui^  =  iyo^/kj^.  The  momentum  of  the 
neutral  atom  beam  exerts  pressure  which  reduces  the  energy  of 
the  waves  in  the  plasma. 

Comparison  of  the  effectiveness  of  the  three  mechanisms  indicates  that 
feedback  gain  required  for  stabilization  by  mechanism  (a),  density  smoothing, 
is  rather  large.  Momentum  control  of  the  neutral  beam,  mechanism  (b),  can  be 
more  effective  than  density  control  if  vx  ~  v£*  At  the  same  time,  mechanism 
(b),  as  equivalent  minimum  B  configuration  may  have  a  relatively  long 


energy  replacement  time  of  the  plasma.  The  alternative  momentum  control, 
mechanism  (c),  can  be  better  than  mechanism  (b)  for  very  small  instability 


growth  rate,  y/uu^,  but  it  is  less  effective  than  mechanism  (a). 

1  2 

In  a  numerical  example,  it  is  assumed  that  a*  =  -  a,.  >  0,  b  =  —  (k.a.) 

G  1  Z  i  1 

=  k2t  —~  =  2(10-2)  or  10-6,  0  =  arg(u;c/uu  )  ,  uu  /uu  =  20  or  200,  v  /r*  =  650 

X  2  ’  °  f  e  s  e  x  e 


MJj 

c 

-  normalize  beam  momentum.  By  introducing  these  values  into  Eqs.  (4.1) 
and  (4.2),  the  calculated  values  of  instability  growth  rate  and  minimum 
feedback  gain  required  for  stabilization  are  obtained.  Figure  11  shows 
the  growth  rate  varying  with  phase  shift  for  various  values  of  the  nor¬ 
malized  feedback  gain  uu^/uu  .  Figure  12  indicates  the  minimum  feedback  gain 

f  e 

I  .  required  for  stabilization  varying  with  the  normalized  neutral 
1  r  c'rnm 

beam  momentum  v  /rju  .  It  is  seen  that  the  value  of  |  i.  / j  I  .  decreases 
x  e  1  f  e 1  m  i  n 

as  .■  />.  increases  when  the  feedback  stabilization  for  resistive  drift 


s  e 


wave  is  concerned. 
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4.3  Modulated  Microwave  Energy  Sources 

The  possibility  of  feedback  stabilization  for  resistive  drift  waves 

of  a  thermonuclear  plasma  by  the  remote  injection  of  neutral  atom  beams 

has  been  analyzed  above.  In  this  section  some  experiments  to  demonstrate 

remote  feedback  control  of  collisional  drift  waves  by  modulated  microwave 

energy  sources  is  discussed.  The  experiment  is  performed  in  a  cesium  plasma 

10  -3 

of  the  Q-machine  with  T  =  0.24  eV,  n  =  5(10  )  cm  and  B  =  4.0  kG. 

e  e 

A  half-wave  dipole  antenna  (located  about  2  cm  outside  the  plasm  olumn) 
is  used  as  a  remote  suppressor,  and  a  Langmuir  probe  is  use  as  ..  detector. 
The  amplified  and  phase-shifted  ion-saturation  current  instability  signal 
is  introduced  to  modulate  the  microwave  energy  source  so  that  the  micro¬ 
wave  power  output  is  proportional  to  the  instability  amplitude.  The 
microwave  energy  heats  the  plasma  by  resonance  absorption.*’*’ 

For  optimum  stabilization,  the  modulated  microwave  power  must  be 
applied  to  the  region  of  maximum  instability  amplitude  and  approximately 
180°  out  of  phase  with  the  local  density  oscillation.  The  measured 
results  from  the  plasma  for  (a)  the  dependence  of  instability  amplitude 
and  frequency  on  feedback  phase  shift  and  gain,  (b)  the  radial  distri¬ 
bution  of  plasma  density,  and  (c)  the  change  of  Magnetic  field  strength 
on  plasma  stabilization  with  microwave  feedback  are  similar  to  those 
obtained  by  the  electrostatic  feedback  (see  Figs.  7-10). 
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4.4  Modulated  Uniformly-Distributed  Electron  Sources 


An  inhomogeneous  plasma  confined  by  a  main  magnetic  field  in  the 
axial  direction  is  considered.  The  particle  density  n  varies  m  the  x- 
direction.  The  feedback  stabilization  with  modulated  electron  sources  uniformly 
distributed  in  the  plasma  is  applied  to  a  very  weakly  collisional  regime 
where  the  electron  mean  free  path  is  greater  than  the  parallel  wavelength. 

In  this  case  the  motion  of  electrons  can  be  described  by  fhe  kinetic  equa¬ 
tion  in  the  drift  approximation.  The  kinetic  equation  for  perturbed  electron 
distributions  with  collision  and  source  terms  in  the  simplest  form  is 


given  by 


where 


if.  E  df  qE  df  df. 

-1£  +  -X  — _  —1  +  v  — i£  =  -„(f  -  n  f)+S 

dt  B  dx  m  dv  z  dz  Ve  ie  le  H  e 

z  z 

(4.3) 


f  =  fff,  \T,  t)  =  the  Boltzmann  distribution  function  (see 
section  2.4) 

f  ,  f.  the  zero-order  and  first-order  distribution 
oe  ie 

functions 

E  =  electric  field  in  the  y-direction 

y 

m,  q  =  the  electron  mass  i  charge 

n  ,  n  ,  v  ,  v  =  the  electron  densities,  velocity,  .and  collision 
oe  ie  e  e  c 

frequency 


f  =  f  /n  =  the  Maxwelliam  velocity  distribution. 
M  oe  oe 


The  feedback  electron  source,  S  ,  is  distributed  throughout  the  plnsra  and 

responds  linearly  to  the  local  electron  density  perturbation.  Based  on 

,  .  .  _  .  .  ,  59,60,69 

this  assumption,  the  source  term  can  be  given  by 


S  =  -  uj  n  f 
e  f  le  M 
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By  the  Fourier  component  analysis  of  Eq.  (4.3)  with  the  aid  of  Eq.  (4.4), 
solving  for  f  and  integrating,  the  ratio  of  the  electron  density  at  local 
perturbed  state  to  the  electron  density  at  initial  equilibrium  state. 


n.  /n  ,  is  obtained 
le  oe 


n. 

ie  _  aJ 


3{ft  1  -  (at  -  -  it/e)  f  [fMdv2/(a)  -  kzvz  4-  iu  ) 

T  1  “  i(^f  -  ^e)  j  ‘  kzVz  +  ivj] 


(4.5) 


in  which 


q0  =  q0  =  work  function,  0  =  scalar  potential, 

T  =  kinetic  temperature,  • 

k  T  . 

.  y _ dn 

®  “  qB  ndx’ 

M  z 

k  ,  k  ,  u>  =  have  been  defined  previously, 
y  z 

In  the  case  of  short  mean  free  path  of  electrons  (for  relatively  strong 

2 

collisions),  in  «  k  v  «  V  and  wv  -  (k  v  )  ,  Eq.  (4.5)  reduces  to 
Z  16  G  G  ZlG 


n-  _  (u*  +  i(k  VT  )  !v. 

ie  _  qje _ z  Te _ ei 

”  T  2 

noe  in  +  ixu  +  i(k  v  )  /y  . 

i  z  Te  ei 


(4.6) 


where  v  =  (2  T/m)2,  electron  thermal  velocity,  and  v  .  =  electron-ion 
Tg  6  x 

collision  frequency.  In  the  case  of  long  mean  free  path  of  electrons 
(for  relatively  weak  collisions),  y  <<r  k  v  and  in  «  k  v  ,  Eq.  (4.5) 

G  ZlG  ZlG 

becomes 


n  M  1  t  i(ir/2)  2  (a)  -  m*  +  iy  )/(|k  Iv  ) 

re  _  gj  _ _ _ _ e  z  1  l  e 

noi  T  1  +  (fl/2)^  (inf  -  ye)/(|kJvTe) 


(4.7) 


If  the  ion  density  is  also  perturbed  by  the  electron  sources,  the  ratio 
of  ion  density  at  local  perturbed  state  to  ion  density  at  initial  equili¬ 
brium  state,  n,,/n  is 

li  oi 

n,  .  ,  , 

1 1  _  -  bm 

n  .  T  in  +  bin 

in 


(4.8) 
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where  b  =  -  k  a.  as  defined  previously. 

For  b  i  tu^/k^v^  «  1,  from  the  quasineutrality  condition  the  dis¬ 
persion  relation  is  obtained. 


(D 


—  U)* 


2b  + 


u) , 


Te, 


+  i(f) 
£ 


UP* 


|  k  I  v. 


Te 


2b 


,1[A  Wf 

2}  ]ir^ 


Te 


(4.9) 


Feedback  stabilization  can  be  achieved  with  purely  imaginary  feed¬ 
back  frequency  uu^  such  that  argluu^/up*)  =  -  90°  (or  +  270°) ,  and  gain 
required  for 


|u)fu)*|  >  2b (id*)2  (4.10) 

Feedback  stabilization  can  also  be  achieved  with  phase  shift  0°  if  the 
condition 


up 


f 


Te 


>  2b 


is  satisfied.  Thereofore,  the  optimum  phase  shift  for  the  feedback  sta¬ 
bilization  will  be  between  -  90°  (or  +  270°)  and  0°.  The  electrostatic 
feedback  by  phase  shift  of  -  90°  and  the  neutral  atom  beam  feedback  by 
phase  shift  of  0°  are  corresponding  to  the  phase  shift  of  the  modulated 
electron  source  feedback.  Although  the  approaches  of  the  feedback  stabili¬ 
zation  are  different,  in  general,  the  feedback  stabilization  by  the  mod¬ 
ulated  electron  sources  is  in  good  agreement  theoretically  and  experi¬ 
mentally  with  the  corresponding  stabilizing  effects  by  the  electrostatic, 
neutral  atom  beam  momentum  or  modulated  microwave  energy  source  feedback. 
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5.  MAIN  COMPARISONS  AND  CONCLUSIONS 


5.1  Main  Comparisons 

In  the  magnetic  field  configuration  stabilization,  the  internal 
axial  magnetic  field,  the  rotational  transform,  and  the  external  axial 
magnetic  field  of  screw  (bumpy  or  diffuse)  pinch  together  with  the 
azimuthal  magnetic  field  to  form  some  fundamental  elements  of  magnetic 
fields  for  plasma  stabilization.  The  magnetic  shear  and  the  magnetic 
well  (including  minimum  B  for  open-end  system  and  minimum  average  B  for 
closed-end  system)  have  shown  good  results  to  control  the  hydromagnetic 
(MHD)  and  drift  instabilities.  The  multipole  configuration  consists  of 
the  linear  and  toroidal  multipoles.  The  linear  multipole  device  has  been 
operated  as  a  quadrupole  in  the  collisionless  regime.  The  toroidal  mul¬ 
tipoles  have  been  studied  from  low  (3  equilibrium  and  flute  instability 
to  high  (3-ballooning  stability.  The  magnetic  shear,  the  magnetic  well, 
and  the  multipoles  can  be  considered  as  the  advanced  magnetic  field 
configuration  for  plasma  stabilization. 

The  high  magnetic  power,  ,  required  per  unit  length  of  the  confinement 
chamber  for  these  configuration  is  approximately  given  by 


P  ~  C 
mag 


B2RT| 

NA 


(5.1) 


where 

B  =  the  magnetic  induction,  or  magnetic  field  strength, 

”  =  the  coil  radius, 

A  =  the  cross-sectional  area  of  the  coil, 

N  =  the  number  of  coil  turns  per  unit  length  of  the  confinement 


chamber , 
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T)  =  the  resistivity  of  the  coils,  and 
C  =  a  constant. 

Among  the  dynamic  stabilization  methods  all  use  oscillating  elec¬ 
tric  field  or  alternating  current,  except  the  electron  temperature  modu¬ 
lation  method,  to  stabilize  the  drift  waves  of  a  thermonuclear  plasma. 

In  the  oscillating  electromagnetic  field  method  and  the  oscillating  non- 
uniform  electric  field  method,  both  electric  and  magnetic  fields  are 
considered.  In  the  electron  temperature  modulation  method, 
microwave  (or  radiofrequency  field)  and  optical  energy  sources  can  be 
employed  to  heat  the  plasma  by  remote  control.  The  oscillating  magnetic 
field  for  use  in  the  dynamic  stabilization  is  usually  produced  by  alter¬ 
nating  current  or  currents. 

Of  the  methods  discussed  above,  *-he  dynamic  stabilization  for  drift 
instability  of  a  collisionless,  inhomogeneous  plasma  is  most  effective. 

The  axial  alternating  current  method  used  for  dynamic  stabilization,  how¬ 
ever,  is  potentially  effective  for  collision-dominated  drift  waves  of 
the  plasma. 

Dynamic  stabilization  has  demonstrated  the  possibility  of  stabilizing 
low-frequency  drift  waves  (or  broad  turbulences)  and  improving  contain¬ 
ment  time  of  a  thermonuclear  plasma.  Apart  from  electron  temperature 
modulation  by  microwave  heating,  modulated  neutral  beam  or  optical  energy 
source  (such  as  weak  laser  beam)  may  be  useful  for  dynamic  stabilization 
of  a  fusion  reactor  by  the  remote  control  and  detection  methods. 

The  power  required  for  dynamic  stabilization,  ,  by  t'10  high- 

,  70 

frequency  electric  field  at  frequency  U  is  given  by 


P,  =  bn  IV,  . 
dyn  o  ■ 
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vhere 

=  equilibrium  plasma  density, 

T,  b  =  as  defined  previously,  and 

bn^T  =  a  fraction  of  the  interrial  (or  thermal)  energy  of 
the  plasma. 

Electrostatic  feedback  for  which  feedback  elements  immersed  in,  or 
in  contact  with,  plasma  have  proved  effective  to  stabilize  plasma  insta¬ 
bilities  and  improve  plasma  confinement.  Particularly,  the  collisional 
drift  instability  in  the  oscillatory  regime  has  been  stabilized  by  immersing 
electrostatic  probes  to  draw  electric  current  from  the  plasma.  Since 
material  probes  cannot  resist  high  temperature  in  a  thermonuclear  plasma, 
the  electrostatic  feedback  method  is  probably  limited  to  laboratory 
plasmas  of  moderate  density  and  temperature. 

The  application  of  feedback  control  to  a  nuclear  fusion  reactor,  how¬ 
ever,  depends  critically  on  the  possibility  of  remote  control  and  detection 
method.  Therefore,  the  remote  control  of  neutral  atom  beam  feedback  to 
stabilize  resistive  drift  waves  and  the  modulated  microwave  energy  sources 
to  stabilize  collisional  drift  waves  has  been  proposed.  Charged  particle' 
beams  and  optical  beams  (including  laser  beam)  can  also  no  used  for  remote' 
control  and  detection  in  the  feedback  stabilization  processes. 

The  effectiveness  of  some  experimental  results  on  the  low-frequency 
drift  waves  and  the  good  qualitative  agreement  with  theoretical  expecta¬ 
tion  indicate  that  the  study  of  feedback  stabilization  (as  well  as  dynamic 
stabilization)  may  be  a  fruitful  approach  to  increase  the  containment 
time  of  a  thermonuclear  plasma. 
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The  power  required  for  linear  feedback  stabilization,  by  the 

parallel  electron  current  is  given  by 


P...  =  a  bn  Tu)r 
fb  o  f 


(5.3) 


in  which  =  uu^  =  instability  frequency,  a  =  q0/T  ~  n/n  =  normalized 
instability  amplitude,  and  bnQT  remains  a  fraction  of  the  internal  (or 
thermal)  energy  of  the  plasma. 

The  ratio  of  power  required  for  the  feedback  stabilization  to  that 
required  for  the  dynamic  stabilization  on  the  basis  of  Eqs.  (5.2)  and 
(5.3)  is 


P 


fb 


/P, 

dyn 


a 


(5.4) 


which  is  much  smaller  than  unity  when  a  «  1  and  uu^  <  0  in  the  most 
practical  cases. 

The  high  magnetic  power  requirement  of  the  advanced  magnetic  lie  Id 
stabilization  system  is  equivalent  to  the  total  magnet  loss  per  unit 
length  of  the  confinement  (or  vacuum)  chamber  due  to  heating,  confinement, 
and  stabilization  of  the  plasma. 

For  purposes  of  comparison  and  conclusions,  the  main  advantages  and 
disadvantages  of  the  plasma  stabilization  systems  are  given  below: 


5.1a  Magnetic  Field  Configuration  Stabilization 

•  Provision  for  plasma  heating,  confinement  and  stabilization 

•  Self-stabilization  in  appropriate  magnetic  field  configuration 

•  Suppression  of  the  hydromagnetic  (or  rillD)  instabilities  in 


particular 
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•  Complex  stabilizing  magnetic  fields  required 

•  Complex  winding  and  high  cost  of  the  solenoidal  conductor 

•  High  power  required  for  magnetic  field  generation  and  operation 

5.1h  Dynamic  Stabilization 

•  Possibility  of  stabilizing  many  modes,  different  low-frequency 
drift  instabilities,  and  broad-based  turbulences  simultaneously 

•  Suppression  of  low-frequency  drift  instabilites  in  particular 

•  Stabilizing  frequency  much  higher  than  instability  frequency 
and  no  phase  shift  needed 

•  Excessive  power  required  to  maintain  new  dynamically  stable 
equilibrium  for  high  growth- rate  instabilities. 

5.1c  Feedback  Stabilization 

•  Relatively  simple  magnetic  geometries  required  when  appropriate 
feedback  and  phase  shift  systems  are  selected. 

•  Lower  power  requirement  for  stabilization  because  when  the 
amplitude  of  instability  is  reduced,  the  feedback  amplitude 
reduces  accordingly. 

•  Feedback  can  interact  only  with  one  unstable  mode  near  marginal 
instability;  it  is  difficult  to  stabilize  many  modes  simultaneously. 

•  Application  of  feedback  control  by  means  of  material  probes  and 
electrodes  limited  only  to  plasmas  of  moderate  density  and 
temperature.  Remote  control  and  detection  methods  are  preferred. 

•  Linear  feedback  (feedback  frequency  instability  frequency 

with  specific  phase  angle)  can  stabilize  different  drift  insta- 
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•  Nonlinear  feedback  (repetition  rate  <  instability  frequency,  with 
specific  phase  angle)  effective  to  suppress  perturbations  in 
the  nonlinear  feedback  constant  corrective  signal  is  applied 
whenever  the  instability  amplitude  grows  above  a  prescribed 
level.  Higher  modes,  however,  are  easily  destabilized. 


5.2  Conclusions 

•  In  the  magnetic  field  configurations,  the  principles  of  magnetic 
shear  and  magnetic  well  will  continue  to  play  an  important  role 
in  plasma  stabilization.  The  trend  of  the  magnetic  shear  and 

the  magnetic  well  is  toward  more  sophisticated  and  more  precisely 
controlled  economic  magnetic  field  geometry  to  stabilize  hydro- 
magnetic  and  drift  instabilities. 

•  Dynamic  stabilization  (0  »  no  phase  requirement)  can 

stabilize  many  modes,  different  low-frequency  drift  instabilities, 
and  broad-band  turbulences  simultaneously .  Excessive  power,  how¬ 
ever,  is  required  to  maintain  new  dynamically  stable  equilibrium 
for  high  growth-rate  instabilities. 

•  Feedback  stabilization  (selected  frequency  and  specific  phase  shift 
required)  can  interact  only  with  a  part  of  the  plasma  and  one 
unstable  mode  near  marginal  instability.  Feedback  control  by 
means  of  material  probes  and  electrodes  is  limited  to  plasmas 

of  moderate  density  and  temperature.  Linear  feedback  (u)f  =  u^) 
can  reduce  drift  instabilities  effectively  and  suppress  the 
m  =  1  mode  of  flute  instability  reasonably.  Nonlinear 
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feedback  (repetition  rate  <  instability  frequency  ug^)  is 
effective  to  suppress  low-frequency  drift  instabilty  at  repetition 
rates  below  instability  frequency  but  higher  modes  are  easily 
destabilized. 

•  There  are  more  basic  differences  between  dynamic  stabilization  and 
feedback  stabilization.  In  the  dynamic  stabilization  the  high 
frequency  (or  radio-frequency)  field  interacts  with  the  entire 
plasma  at  a  new  equilibrium  density  nQ  while  in  the  feedback 
stabilization  the  plasma  equilibrium  remains  unchanged  and  the 
time-varying  suppression  signal  interacts  only  with  the  unstable 
part  of  the  plasma.  Both  dynamic  stabilization  and  feedback  sta¬ 
bilization,  however,  will  prefer  to  have  remote  control  and 
detection  methods  for  use  in  a  controlled  fusion  reactor,  even  a 
rather  advanced  and  costly  technique  is  required  for  remote 
dynamic  or  feedback  control. 

•  At  present,  plasma  self-stabilization  in  appropriate  economic 
magnetic  field  configurations  is  still  uncertain.  If  the  plasma 
self-stabilization  can  prove  to  be  fully  effective  in  the  labora¬ 
tory  scale,  the  actual  design  of  a  controlled  fusion  reactor  with 
the  provision  of  complex  stabilizing  magnetic  fields  will  re¬ 
quire  some  reliable  scaling  laws  to  determine  the  reactor  size. 

The  reliable  scaling  laws  for  controlled  fusion  reactor  design, 
however,  have  not  satisfactorily  been  found.  Since  the  magnetic 
field  configuration  stabilization  is  mostly  effective  in  sta¬ 
bilizing  the  hydromagnetic  instability.  The  dynamic  or  feed¬ 
back  stabilization  which  can  stabilize  the  drift  (.gross  or 
universal)  instability  of  a  thermonuclear  plasma  may  contribute 


to  controlled  fusion  research. 
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•  It  is  expected  that  the  research  and  development  of  remote  con¬ 
trol  and  detection  method  for  potential  dynamic  stabilization 
or  feedback  stabilization  of  a  thermonuclear  plasma  will  pro¬ 
bably  become  important  in  developing  a  controlled  fusion  reactor. 
Either  the  dynamic  stabilization  or  the  feedback  stabilization 
may  be  combined  with  the  magnetic  field  configuration  to  sta¬ 
bilize  a  thermonuclear  plasma  so  that  the  minimum  requirement 
for  a  controlled  fusion  reactor  may  be  attained  and  satisfied. 

In  the  preceding  discussion  and  analysis  a  number  of  new  ideas  which 

can  be  developed  into  original  methods  for  understanding  plasma  instability 

and  improving  plasma  stabilization  have  aeen  introduced.  The  main  fields 

of  the  new  ideas  may  consist  of  the  combination  of  magnetic  shear  and 

magnetic  well,  and  the  combination  of  magnetic  field  configuration  and 
dynamic  or  feedback  stabilization  with  remote  control  and  detection  methods. 

In  the  main  fields,  some  specific  ideas  such  as  the  combined  effect  of 
rotational  transform, magnetic  shear  and  minimum  B  on  a  toroidal  plasma, 
extention  of  the  present  magnetic  shear  analysis,  dynamic  stabilization 
by  remote  microwave  (cr  radiofrequency)  field  control  (see  Fig.  6),  etc., 
can  be  developed  to  make  contributions  for  a  controlled  fusion  reactor. 
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